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ABSTRACT 

Observations have established that the diffuse intergalactic medium (IGM) at z ~ 3 is emiched to ~ 
10~^'^ solar metallicity and that the hot gas in large clusters of galaxies (IGM) is enriched to 1/3 — 1/2 Zq 
at z = 0. Metals in the IGM may have been removed from galaxies (in which they presumably form) 
during dynamical encounters between galaxies, by ram-pressure stripping, by supernova-driven winds, 
or as radiation-pressure driven dust efflux. This study develops a method of investigating the chemical 
enrichment of the IGM and of galaxies, using already completed cosmological simulations. To these 
simulations, we add dust and (gaseous) metals assuming instantaneous recycling, and distributing the 
dust and metals in the gas according to three simple parameterized prescriptions, one for each enrichment 
mechanism. These prescriptions are formulated to capture the basic ejection physics, and calibrated when 
possible with empirical data. Our method allows exploration of a large number of models, yet for each 
model yields a specific (not statistical) realization of the cosmic metal distribution that can be compared 
in detail to observations. Our results indicate that dynamical removal of metals from ^ 10^'^ M0 
galaxies cannot account for the observed metallicity of low-column density Lya absorbers, and that 
dynamical removal from ^ lO^^'^M© galaxies cannot account for the IGM metallicities. Dynamical 
removal also fails to produce a strong enough mass- metallicity relation in galaxies. In contrast, either 
wind or radiation-pressure ejection of metals from relatively large galaxies can plausibly account for all 
three sets of observations (though it is unclear whether metals can be distributed uniformly enough in 
the low-density regions without overly disturbing the IGM, and whether clusters can be enriched quite 
as much as observed). We investigate in detail how our results change with variations in our assumed 
parameters, and how results for the different ejection processes compare. 

Subject headings: cosmology: theory - intergalactic medium - galaxies: formation, starburst - 
hydrodynamics 



1. INTRODUCTION 

In the standard hot big-bang cosmological model, ele- 
ments of atomic number Z > 2 cannot form in appreciable 
quantities until the first stars form. Thereafter, the uni- 
verse becomes progressively enriched with heavy elements 
(metals) as stars release these fusion products in stellar 
winds or supernova ejecta. Stars only form efficiently in 
dense gas, almost all of which is bound in galaxies (or pro- 
togalaxies at high z) with deep gravitational wells. Despite 
this, the observed intergalactic medium (IGM) shows sub- 
stantial metal enrichment at all redshifts and at all densi- 
ties yet measured, from the lowest-density Lya absorbers 
at 2; > 3 with ^ 1/1000 solar metallicity, to the hot 
gas with 1/3-1/2 Zq bound in present-day galaxy clus- 
ters. A significant fraction, perhaps even a majority, of 
cosmic metal appears to lie in the IGM. At very high red- 
shift {z 3> 5) population III stars could enrich the IGM 
to a low level (e.g., Carr, Bond & Arnett 1984; Ostriker 
& Gnedin 1996; Haiman & Loeb 1997; Tegmark et al. 
1997; Abel et al. 1998), but since most cosmic metals 



presumably form in stars in galaxies, a fundamental ques- 
tion arises as to how these metals escape their progenitor 
galaxies and spread throughout the IGM. 

Efforts to understand the enrichment of the IGM by 
galactic stars (after any Pop. HI epoch) have focused on 
three mechanisms whereby metals could be removed from 
a galaxy. First, metal-enriched gas (or stars that later ex- 
plode as Type la supernovae) might be unbound during 
a merger or tidal interaction with another galaxy, or by 
the ram pressure of the IGM through which the galaxy 
moves. We shall combine these processes under the name 
of 'dynamical removal'. Second, the energy input from su- 
pernovae may impart sufficient kinetic and thermal energy 
to galactic gas for it to escape the gravitational well of the 
galaxy. We shall denote this process as the 'galactic wind' 
mechanism. Third, the radiation pressure on dust grains 
due to stellar light may exceed the gravitational force of 
the matter, leading (if the dust can decouple from the 
gas) to an outflow of dust. We will denote this possibility 
by 'radiation-pressure ejection' or 'dust ejection'. -'^ (The 
ejected dust adds metals to the intergalactic gas when it 



^Dust will, however, also be ejected in dynamical removal or wind ejection along with the gas. 
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is destroyed by thermal sputtering.) 

Investigation of the metal enrichment of the IGM re- 
quires not only some understanding of how these mecha- 
nisms function in a given galaxy, but also calls for knowl- 
edge of the properties, distribution, and evolution of galax- 
ies in a cosmological context. Studies of these matters have 
generally adopted one of two rather different, yet com- 
plementary, approaches: the 'numerical' and the 'semi- 
analytic' methods. 

In the numerical approach, an attempt is made to nu- 
merically evolve the state of a single sample of the universe 
from some set of initial conditions, by solving equations 
of motion encapsulating the most relevant (gravitational, 
hydrodynamical, etc.) physics. Among the advantages 
of this approach are that it captures the incorporated 
physics very well, and that it generates a specific real- 
ization of a possible cosmological volume. On the other 
hand, the available computing power limits the amount 
of detailed physics and the dynamic range that can be 
handled. Also, some physical processes must be 'param- 
eterized', and the large computational time prohibits the 
investigation of large regions of parameter space. The nu- 
merical method has been used by Gnedin (1998) to study 
the enrichment of the IGM by winds and dynamics and 
employed by Gen & Ostriker (1999) and Gnedin & Ostriker 
(1997) to study IGM enrichment by dynamics or other (un- 
specified) processes. Cluster enrichment by galactic winds 
has been studied using simulations by Metzler & Evrard 
(1994; 1997) and Murakami & Babul (1999), and the en- 
richment of the intracluster medium (IGM) through dy- 
namics has been numerically studied by Abadi, Moore & 
Bower (1999), by Balsara, Livio & O'Dea (1994), and by 
Quihs, Moore & Bower (2000). 

In the semi-analytic approach, predictions are made on 
a statistical basis, by layering together a number of pre- 
scriptions that are individually derived either from theory 
or observation. This approach allows the investigation of 
large regions of parameter space, and employs somewhat 
more complicated physical prescriptions than the numer- 
ical method. On the other hand, the flexibility in the 
input parameters and physical prescriptions leads to a 
corresponding range in actual predictions, and it is not 
always clear which physical processes are accurately cap- 
tured and which are not; nor is it clear that the param- 
eters chosen to best fit the observations are unique. The 
semi-analytic method has been used by a number of in- 
vestigators to address IGM enrichment in greatly varying 
levels of complexity. Nath & Trentham (1997), Ferrara, 
Pettini & Shchekinov (2000), and Madau, Ferrara & Rees 
(2000) have studied IGM enrichment by winds in this way, 
attempting also to calculate the statistical properties of 
the metal distribution in the IGM. Cluster enrichment by 
winds has been studied by Dekel & Silk (1986), Nath & 
Chiba (1995), and David, Forman & Jones (1990), among 
others. The dynamical enrichment of clusters has been 
similarly examined a number of times (e.g., Gunn & Gott 
1972; Renzini et al. 1993). In both cases the overall de- 
gree of enrichment has been assessed, but not the distri- 
bution of the metals. Finally, ejection of dust by radiation 
pressure has been calculated for individual sample galax- 
ies (Chiao & Wickramasinghe 1972; Ferrara et al. 1990; 
Shustov & Vibe 1995; Davies et al. 1998; Simonsen & 
Hannestad 1999), but only rough estimates of the overall 



ensuing IGM enrichment are offered. 

The method used in the present work combines as- 
pects of both approaches. The time-dependent distribu- 
tion of dark matter, gas and stars is taken from already 
completed cosmological smoothed-particle hydrodynamics 
(SPH) simulations computed using the method described 
by Katz, Weinberg & Hernquist (1996). To the gas and 
star particles, metals are added using fairly simple pre- 
scriptions formulated to capture the basic physics of the 
various processes that can transport metals from galax- 
ies into the IGM. This method yields a numerical real- 
ization of the actual distribution of metals in the IGM 
(as well as in galaxies), yet uses relatively little comput- 
ing time so that a large number of models can be tested, 
and the effects of changing both physical prescriptions and 
parameters can be explored. This method has a number 
of limitations that will be discussed in detail in the fol- 
lowing sections, but it can nevertheless yield important 
insight into metal ejection and distribution unobtainable 
by either 'purely' numerical computations (using present 
technology), or semi-analytical methods. 

This paper presents an explication and investigation of 
the method we have developed. We do not attempt to 'fit' 
a set of observations using our results, but rather make 
specific assumptions about our input parameters based 
on independent considerations, and compare the predic- 
tions to observations where available. We also investigate 
directly the variations in predictions resulting from differ- 
ent parameter choices, and address various methodological 
and numerical considerations. Studies addressing specific 
astrophysical questions using the method presented here 
will be pubhshed separately (e.g., Aguirre et al. 2001a, b). 

We have organized this paper as follows. Section || de- 
scribes our calculation method, with subsections |2.l|, b.2L 



and 2.3 detailing prescription for dynamical, wind, and ra- 
diation pressure ejection of metals, respectively. These sec- 
tions also describe the parameters used in the calculation 
and the fiducial values for these parameters (summarized 
in Tables |l], ^, || and |^) . In § || we review observations of 
cluster, galaxy, and Lya forest metallicity with which we 
will compare our results. Section |^ briefly compares the 
SPH simulations (to which we a ppl y our method in this 
paper) to observations. Sections ^3|, ^ and ^ present and 
discuss the trials we have run and the results obtained. 
We summarize our findings and draw general conclusions 



2. METHOD 

The procedure begins with a set of coarsely-spaced (ev- 
ery 10^ — 10^ yr) snapshots from an SPH cosmological sim- 
ulation, each containing the states of the dark matter, star, 
and gas particles at a given redshift. We shall index these 
particles by 'i', and superscript them by particle type as 
'd', 's' or 'g' respectively. We will also denote by tc the 
time of the output being processed, and by tp the time at 
the previous output. 

Beginning with the first output time where stars exist, 
the procedure for each time step is as follows. For a star 
particle i, the stellar mass created since the last step is 
Amf{tc) = ■mf{tc) — mf{tp). This mass was taken from a 
gas particle also indexed by «, so the first task is to trans- 
fer the metals associated with this gas mass from the gas 
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particle to the star particle: 



An^it,) X K(tp)/mf(tp)] 
Amtit,) X K(ip)/mf(tp)], 



(1) 



the metal mass within a shell of width dr at radius r is 
proportional to W^'^^{r,h'^™'^{9,(j)))dr; uniform distribu- 
tion of metals within the angular region thus corresponds 
to M^g'P(r,/i™''^(6',0)) 
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where Wi is the metal mass of the zth particle. 

We next assume that each unit of forming stellar mass 
instantaneously ejects units of metal mass into the gas, 
so that the star particle i adds a metal mass of Am* {tc)xy^, 
to the gas particles. The metals are distributed over the 
gas particles in a manner appropriate to the particular pro- 
cess, as described in the next few sections. After all metal 
is deposited, the process repeats for the next simulation 
output time. 

2.1. Prescription 1: Local Metal Distribution 

The first prescription for metal distribution is the sim- 
plest. The metal mass Amf (tc) x is 'scattered' over the 
32 gas particles nearest to the star particle i, weighted by 
the SPH smoothing kernel W{r, h), where 2h is the radius 
of the sphere (about particle i) containing exactly 32 gas 
particles, and r is the distance from particle i to the gas 
particle in question (see Hernquist & Katz 1989). Since 
stars form in these simulations only in dense, cool regions 
with a converging flow, this prescription places metals only 
in bound regions such as galaxies.^ Thus with this pre- 
scription intergalactic metal must leave galaxies by being 
carried by a gas particle that is dynamically removed from 
a bound region. 

2.2. Prescription 2: Ejection by Galactic Winds 

Galactic winds can distribute metals 'non- locally', i.e. 
disperse them into gas far away from where they are 
formed. The prescription used to model this physical ef- 
fect is as follows. First, we divide the gas and star parti- 
cles into bound groups (i.e. galaxies) using the SKID algo- 
rithm.^ At each simulation step, each bound group is con- 
sidered in turn. As described above, the new metal mass 
Amf{tc) X in each star particle i within the group is 
computed. A fraction (1 — Ycj) of this metal is distributed 
among the 32 gas particles nearest to star particle i as 
per prescription 1. The remaining fraction is added to the 
tally for the new metal in the group, Aw^'^^: 



for each galaxy requires some understanding of the physics 
of galactic winds. The idea that galaxies might drive 
outflowing winds has a fairly long history. Galactic-scale 
winds, in which supernova bubbles can overlap and drive a 
coherent wind across the galaxy before they can cool, have 
been theoretically investigated in dwarf (e.g., Mac Low 
& Ferrara 1999), elliptical (e.g., David, Forman & Jones 
1990; 1991 and references therein), and starburst spiral 
(e.g., Heckman et al. 2000, hereafter HLSA; Lehnert & 
Heckman 1996) galaxies. In the last case, such winds can 
be observed in some detail as high velocity bipolar outflows 
seen in many starburst galaxies such as M82 (e.g., Lehnert, 
Heckman & Weaver 1999). Even if the conditions to drive 
a coherent galactic-scale wind do not exist, winds capable 
of ejecting matter into the IGM may still develop: a two- 
phase ISM with a hot phase fed by supernova remnants can 
lead to a stochastic 'steady-state' in which some fraction 
of matter has high enough kinetic energy - perhaps aided 
by cosmic-ray pressure (e.g., Breitschwerdt et al. 1991) - 
to escape the galaxy (especially near concentrations of su- 
pernovae). Winds of this sort have been investigated the- 
oretically by Efstathiou (2000), Ferrara & Tolstoy (2000), 
Breitschwerdt & Schmutzler (1999), and Ferrara, Pettini 
& Shchekinov (2000). Thus we see that winds may be 
starburst-driven, driven by 'quiescent' star formation, or 
cosmic ray-driven, and that they may be 'global' or 'lo- 
cal'. These distinctions can readily break down, however: 
in dwarfs or in galaxies with very rapid star-formation, the 
local/global distinction breaks down, and at high- 2: (where 
star formation is vigorous and mergers common) there is 
probably no clear line between starbursts and quiescent 
star formation. 

Despite their differences, all of the wind types we have 
described share some common features: 
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(2) 



1. The energy release in supernovae is the ultimate 
source of the wind energy. Some critical supernova 
rate is necessary for a galactic-scale wind to form 
and blow out of the disk. 

2. The wind speed may exceed the escape velocity of 
the progenitor galaxy, but the wind and the swept- 
up material must stall or become pressure-confined 
at some radius. 

3. It is physically reasonable for the wind's energy to 
be tied to the star formation rate. 

The method used in our investigation is based primarily 
on observations of galactic-scale 'superwinds' and has been 
formulated to capture these three key physical features in 
a simple and general way, so that it can be reasonably ap- 
plied to galactic winds of all types (with varying degrees 
of confidence). First, we assume some critical SFR/ (area), 
SFRcritj below which wind development is suppressed. For 

■^If stars form in bound groups with < 32 gas particles, some metals will be placed in particles outside the bound group. 

^SKID is publicly available at ittp://www-hpcc. astro. Washington. edu/tools . We require also that these groups have at least 4 (gas+star) 



We then distribute the metal mass w^'p within a radius 
^wmdj-^^^-j g^ijout the 'center of star formation' 



(3) 



where xf is the position of star i. This simulates the non- 
local dispersal of metal into the regions where winds carry 
and deposit them. The metal is distributed so that within 
the angular ranges [cos 61, cos 6* + Acose] and [0,0 + A^], 
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particles and a minimal overdensity ot 5U; groups ol between 1 and 4 particles are treated as 'ungrouped', but essentially no star formation 
takes place in these areas. 
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thermal winds, this should physically correspond to a rate 
above which supernova remnants can overlap before cool- 
ing (David, Forman & Jones 1990; Efstathiou 2000; Heck- 
man, Armus & Miley 1990). This assumption is supported 
by observations indicating that superwinds in spirals de- 
velop when the SFR/unit area (averaged over the disk) 
reaches a critical value (Martin 1999; Heckman 2000). 
Since the 'areas' of spiral galaxies are not robustly deter- 
mined in the simulations, we compute galaxy areas from 
their masses, using the empirical mass-radius relation^ of 
Gavazzi, Pierini & BoseUi (1996), who find (area)oc M^-'^^. 

We allow all galaxies to drive winds, with an initial wind 
velocity Wout and mass outflow rate mout- Guided by ob- 
servations of starburst-driven superwinds (Heckman et al. 
2000; Martin 1999) indicating mass outflow rates similar 
to the galaxies' SFRs (for a wide range of SFRs), we ex- 
press rhout in units of the SFR.^ Assuming that one su- 
pernova forms per AfioolOO Mq of star formation and re- 
leases 10'^^i?5i ergs of energy, the outflow rate m x (SFR) 
(in M0/yr) can be related to the fraction x of the super- 
nova energy that is incorporated into the wind's kinetic 
energy mvl^^ x (SFR), by 



X 



E^i Afioo, 



(4) 



where Wout is the wind velocity. 

We use fixed values of SFRdit and x that are based 
on the data available for superwinds (see 



2.2.1 



be- 
low). For Wout, we assume a uniform distribution of val- 
ues with mean v^^^ and width CTout- Given x and Vont 
for a given galaxy, we compute rhout using Eq. |^. For 
galaxies with SFR/ (area) < SFRdit, we attenuate Vont 
by a factor (SFR/SFRcrit)'' and attenuate x by a fac- 
tor (SFR/SFRcrit)^''- This lowers both the energy of 
the wind and the energy of the initial shell by a factor 
(SFR/SFRcrit)'^'', and therefore suppresses winds in low- 
SFR galaxies as desired. 

We assume that in each direction {9, (j>) the wind flows 
to some maximum radius /i"'™'^(0, </>), then joins the dy- 
namics of the ambient gas. This radius, within which the 
galaxy's metals are distributed, is again determined from 
the physics of winds. When a bubble of supernova-heated 
gas forms and begins to expand, several things can stop 
its growth. First, the ISM of the host galaxy will be swept 
into a shell which, if massive enough, can conflne the wind 
to the galaxy, as happens in an ordinary single supernova 
explosion. If the bubble has enough energy, the shell can 
be blown out of the disk, whereupon it may partially frag- 
ment due to the Raleigh- Taylor instability or because of 
density inhomogeneities in the ambient medium; the hot 
wind can then stream past the shell, entraining it (and 
perhaps other portions of the ISM) into a mass-loaded 
outflow. The process may then repeat. 

For each of Na directions {0,(f>), we compute the maxi- 



mum radius to which the wind could expand by following 
the dynamics of a 'test' shell of physical radius r about 
Tc, with mass m(r) and outward radial velocity v{r). This 
shell represents the initial wave of matter swept up by 
the developing wind, and feels four forces: the wind's ram 
pressure, gravity, the thermal pressure of the ambient gas, 
and the ram pressure of any infalling ambient gas^ Under 
these four forces the shell momentum evolves according to 
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and the mass evolves as 
d 
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d(t> 
dr 
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r^P (Wrad + Hr) (v - Urad - Hr), 



dt 



(m) 



1 



An 



r^p{v - Wiad - Hr) . 

(6) 

Here, p{r), p(r), (f>{r) and Wrad('') are the average density, 
thermal pressure, gravitational potential and outward ra- 
dial peculiar velocity of the ambient medium at radius r, 
Vu}{t) is the wind velocity, and H{z) is the Hubble constant 
at redshift z: 

H^{z) = Hi [(1 + zfVl + f^A + (1 + zf{l -n- ^a)] . 

The averages are done over all particles within the angular 
ranges [cos 9, cos 9 + Acos e] and [0, (j) + A^] and the radial 
range'' [0.95r, 1.05r]. The angles are spaced in Ng seg- 
ments of Acose = '^/Ng in cos 6* and N^ — 2Ng segments 
of A^ = 2tt/N^ in 0, giving Na — NgN^ portions of equal 
solid angle. We choose Na for each galaxy to be l/16th of 
the number of gas-|-star particles in the galaxy; tests show 
that this gives a small enough number of angles for the 
radial integrations to be accurate (see also § ||). 

The entrainment parameter eent is the fraction of am- 
bient material we assume to be swept up by the wind; 
a perfectly homogeneous shell expanding into a homoge- 
neous medium would give eent — 1 , whereas Cent ^ 1 might 
describe a very clumpy ambient medium with a very small 
filling factor that the wind can easily stream past (while 
filling in the 'holes' after passing), or a shell that repeat- 
edly fragments, leaving behind a fraction (1 — eent) of its 
mass. (Note that we have not included e in the terms ac- 
counting for the wind on the shell, as would be appropriate 
if the shell itself had a small filling factor.) 

The wind velocity (r) is Uout , attenuated by gravity: 



,{r) = J {v, 



t)^ 



2[<^(r)-0(ro)], 



(7) 



where tq is the initial radius of the test shell, defined be- 
low. 

The physical radial velocity of the shell with respect to 
the ambient gas is w — Urad ~ Hr, so evolving (mv) and 

*The masses derived by Gavazzi et al. from rotation curves will include some dark matter contribution not inr liiHed in the corresponding 



2.2 



simulation galaxy mass. This should not affect large galaxies too much but may be important in dwarfs. See 

^This mass outflow rate, using the current methodology, does not really determine how much metal leaves the galaxy, merely the 'strength' 
of the wind. 

^This purely momentum-based method is similar to that of Theuns, Mo, & Schaye 2001; other studies such as Tegmark, Silk, & Evrard 
1993, Nath & Trentham (1997) and Scannepieco & Broadhurst 2001 use slightly different approach in which the shell is driven by thermal (not 
kinetic) energy so that the internal energy of the 'bubble' must be evolved. 

'^Other prescriptions for the radius to average over could be used, but experimentation shows that the integration is not sensitive to the 
averaging scheme so we employ this simple one. 
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m until v{r) — Uiad — Hr < 5 x Hr [6 <C 1) gives the 
radius rgtaii and time Tstaii &t which the sheU staUs and 
mixes with the ambient gas (i.e. moves a distance 5 x r 
in another Hubble time; we choose 5 — 0.05). We set 
initial conditions for m, v and r by choosing a radius 
to include a fixed fraction ^ of the galaxy's mass Mgai. 
The initial shell mass is then r7i(ro) — Ccnt^-^gaiZ-^aj and 
viro) = Wquj. (The fraction ^ is calibrated to give values 
of ro similar to the radii a t whi ch winds from starburst 
galaxies are observed; see § |2.2.l| below) . 

The stalling radius must be modified by two effects re- 
sulting from our general method, in which we deposit 
metal generated between times tp and in the IGM 
at time tc rather than at the more appropriate average 
time (tp + tc)/2 + Tstaii- First, if (tp + tc)/2 + Tstaii ex- 
ceeds the time iobs at which the results for metal en- 
richment are desired, we would allow metals to move to 
erroneously large radii; thus we limit the shell to have 
Tstaii < tobs-(ip + ic)/2. Second, if Tstaii + (tp + tc)/2 > 
then metal deposited at will effectively be artificially 
carried by the movement of the gas particles between tc 
and Tstaii + {tp + tc)/2. For example, if the shell stalls in 
a region with radial velocity Wrad and its metal were de- 
posited at tc, then at the appropriate distribution time 
''"stall + (^p + tc)/2 the metal would be in particles at 
radius (roughly) retail + frad [tstaii + {tp - ic)/2], rather 
than at Tstaii- To compensate for this, after the shell 
reaches rgtaii we continue to integrate in radius, but us- 
ing dr/dt = -Vre^dir) for a time Tstaii + {tp - tc)/2 (or 
^obs ~ tc if this is smaller), to reach a final radius ft-wind- 
Then /iwind will be such that metal deposited at /iwind 
will end up at rgtaii after being carried by the gas parti- 
cle movement for time Tstaii + (^p — tc)/2. This assumes 
that the radial peculiar velocity field in the galaxy's neigh- 
borhood changes on a timescale somewhat longer than 
Tstaii + (^p — tc)/2. A similar (small) error of the oppo- 
site sign arises if Tstaii + {tp + tc)/2 < tc- In this case, we 
compensate for this error by integrating dr/dt — Vra,d{r) 
after stalling, for the time {tc — tp)/2 — Tgtaii to obtain 

^wind ■ 

The metal mass is then distributed within the radius 
= niin(rstaii,w), if ^staii > 2ro. If r'''^^" < 2ro, 
the wind has been efficiently confined, and we instead dis- 
tribute the metal 'locally' as per prescription 1. 

The outlined prescription is based on empirical data 
concerning galactic-scale superwinds driven by starburst 
nuclei. For other types of winds (cosmic ray driven winds, 
'quiescent' winds, 'local' winds, etc.) there is far less use- 
ful observational data; most of our understanding of such 
winds derives from theoretical work (e.g., Breitschwerdt 
et al. 1991; Efstathiou 2000, Ferrara & Tolstoy 2000; 
Breitschwerdt & Schmutzler 1999, and Ferrara, Pettini & 
Shchekinov 2000). Rather than invent a new prescription 
for these winds (based on the theoretical work), we have 
chosen to model them with the same prescription as for 
starburst-driven superwinds, but with different parame- 
ters. We shall lump these other types of winds into the 
category of 'quiescent winds', which will be characterized 
by a lower critical SFR (so that essentially all galaxies 
drive winds) , but also a lower x since the winds would not 
be as efficient. 

The general wind prescription captures the essential 
physics of winds well, and and should provide a good pre- 



diction of the dispersal of metals by winds with a given as- 
sumed set of physical parameters. Since we do not model 
the winds in full physical detail, however, a number of 
methodological details merit attention. 

1. Adding metals after the completion of the simula- 
tion precludes any self-consistent treatment of the 
effects of the winds on the gas near the galaxy. This 
means that we cannot assess the importance of ear- 
lier outfiows in the escape of later outfiows and moti- 
vates our approach of assuming a 'steady state' wind 
running into an undisturbed IGM. Also, we do not 
treat the interaction of winds from a given galaxy 
with those from nearby galaxies. Thus galaxies can 
(unrealistically) pollute wind-driving neighbors, but 
only at a small level. 

2. The method also glosses over the detailed structure 
of the outflowing wind; we have lumped the details 
of the wind's interactions with the ambient medium 
into the parameter ecnt7 assumed not to vary with 
radius or among different galaxies. 

3. We neglect the effect of the matter within the test 
shell radius on the wind that drives the shell. This 
is correct for Ccnt = 1 (since then there is no mat- 
ter within the shell) and for Cent = (since then 
the wind passes through the matter). For inter- 
mediate values, the wind will have a somewhat 
lower mass to and higher velocity Vw than the mass 
and velocity m{ecnt),Vw{^cnt) that it should have 
had. Fortunately these effects largely cancel out 
since (in the limit where gravity is unimportant), 
iTT'{ient)vw{<^cnt) — mv^ by momentum conservation. 

4. We assume a purely kinetic wind impinging upon 
the shell at all stages. This leads to the neglect of 
the pressure inside the shell in Equation |^. It also 
allows us to disregard shocks and cooling, since the 
internal energy of the shell does not affect its propa- 
gation. We retain the external pressure term so that 
environments of very high pressure (i.e. clusters) are 
not treated incorrectly. 

5. Although we compute the maximum wind radius 
with good accuracy, the deposition profile within this 
radius, W^^'^{r,h), must simply be assumed. (But 
we will demonstrate that the details of the assumed 
form are not too important). 

6. Our method assumes that the properties of galax- 
ies driving winds changes on a timescale that is long 
compared to the outfiow time. In particular, we as- 
sume that the SFR is constant while the wind prop- 
agates. We check this assumption in § 

7. The SFRs we use are averaged over the interval 
between simulation output times (typically 10* — 
10^ yr). Weinberg et al. (2000) find that such aver- 
ages are a good approximation to the computed 'in- 
stantaneous' SFRs. However, the simulations cannot 
resolve the small (few hundred pc) scale of starburst 
nuclei, so real starbursts may be shorter and more 
intense than the smoothly-varying SFR the simula- 
tions would suggest, and lead to a larger scatter in 
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the SFR/area than occurs in the simulated galaxies 
(see Dave et al. 1999b). This is not a problem for 
more 'quiescent' winds but complicates the use and 
choice of a critical wind-driving SFR if winds are 
assumed to be primarily starburst-driven. 

8. Our method concentrates all star formation at a 
point, from which the wind emanates. This max- 
imizes the effect of the wind and is appropriate 
in modeling starbursts, but overestimates the wind 
effectiveness if star formation is distributed over 
the galaxy. On the other hand, we assume that 
the wind energy is distributed isotropically. Real- 
istically, winds in their early thermally-dominated 
phase will be funneled into a bi-conical outflow (e.g., 
Mac Low & Ferrara 1999) observations show bi- 
conical outflows with solid angles of ^ O.Stt — 2.47r 
radians at small (few kpc) radii (HLSA; Lchnert & 
Heckman 1996). Thus realistic winds would have a 
higher energy per unit solid angle and a more effec- 
tive outflow in those directions. 

9. If ft,*'""^ is small compared to the size of the galaxy 
but Tstaii > 2ro, the given prescription will artifi- 
cially concentrate metals in the radius /i^'"'^ about 
Tc- This indicates that metallicity gradients in the 
objects should not be trusted, and that the efficiency 
of dynamical removal (which cannot easily remove 
metals from galaxy cores) may be suppressed. 

2.2.1. Fiducial Parameter Values for Wind Ejection 

The key parameters we need to determine are the mean 
outflow velocity v^'^^, the width (Tout of the velocity distri- 
bution, the fraction of supernova energy in the wind X; the 
critical SFR/kpc^ (SFRcrit), the enclosed mass fraction ^ 
determining the initial shell radius, and the entrainment 
fraction font- 

We base these values as much as possible on observa- 
tions of galaxies with supernova-driven winds, as recently 
compfled by HLSA and by Martin (1999). HLSA mea- 
sure both the width W ~ 300 — 600 km s^^ and velocity 
offset Vont ~ 100 — 300kms~^ (from the galaxy's inferred 
center-of-mass velocity) of Nal absorption lines. Their in- 
terpretation of these values is that cool material is being 
accelerated by a hot outflowing wind. Vont is the character- 
istic velocity of this outflowing material at small radii and 
iJterm = 'i^out + W/2 ~ 600kms~^ is the inferred 'terminal 
velocity' to which the dense gas is accelerated, i.e. the ve- 
locity where it is roughly comoving with the hot gas (see 
HLSA). This view is supported by the rough agreement 
of the inferred wtorm values with those inferred from X-ray 
data (HLSA; Martin 1999). ^From the combined data set 
of HLSA we take = 600kms-\ tJout = 200 kms"^ 
Very similar values of t^out and W are found for Lyman- 
break galaxies at z 3 (Pettini et al. 2000). Somewhat 
unexpectedly, according to current data the outflow veloc- 
ities do not seem to be correlated with either the SFR or 
the mass of the host galaxies (HLSA; Martin 1999; Heck- 
man, private communication). Hence we use a fixed value, 
though a different prescription may be called for if future 
observations reveal some dependence on galaxy properties. 

^Indeed, as shown in Fig. ^ the wind velocity substantially exceeds 
^Martin (1999) finds a much smaller threshold of (few) X 10"" Mq y 



In using the observed outflow velocities, we are making 
an important assumption: a single velocity characterizes 
the outflowing cool gas and the hotter wind. It is likely, 
however, that the multiple components of the outflow have 
rather different velocities.* Two-component plasma fits to 
the X-ray data tend to give temperatures corresponding to 
velocities that bracket the Wtcrm inferred from absorption 
lines, but the X-ray observations only measure the ther- 
mal energy of the gas; numerical simulations suggest that 
it may have a kinetic energy ~ 2 — 3 times higher, i.e. ve- 
locities of ~ 800- 1000 kms-i (HLSA; see also Strickland 
& Stevens 2000) . If this hotter gas contains a mass compa- 
rable to the cool material, it may dominate the pressure of 
the outflowing gas. Another possibility is that the outflow- 
ing cool gas may be better characterized by the observed 
velocity offset Wout ~ 200—300 kms~^. This would follow if 
the absorption occurred in a thin shell that might break up 
but is inefficiently accelerated by the hotter wind. In this 
case, the ram pressure of the outflow might be dominated 
by the high- velocity but low-density wind, or alternatively 
by the higher-density but lower-velocity entrained clumps. 
This ambiguity is a significant source of uncertainty, but 
the results should be bracketed by models with velocities 
of 300kms~^ and ~ 1000 km s~^, and we will test such 
values in our calculations. 

The mass outflow rate from a wind-producing galaxy 
can be roughly estimated by measuring the column den- 
sity of the wind material, then assuming either a thin shell 
(Pettini et al. 2000; Martin 1999) or a spherical mass- 
conserving wind (e.g., HLSA). Either way, estimated mass 
outflow rates are ^1 — 4 times the galaxy's estimated 
SFR. Using Equation |4[ this can be converted into a fldu- 
cial value of x ~ 1- (This does not mean that we assume 
supernovae drive winds with perfect efficiency, since prob- 
ably E51 7^ 1 or Mioo ^ 1 in Eq. ^ x is proportional to the 
true wind-driving efficiency, calibrated by the observed ve- 
locities and mass outflow rates for winds.) As noted above, 
for SFR/ (area) < SFRcrit we attenuate the wind energy 
by (SFR/SFRcrit)^'^. We use /? = 2 as a default, but in- 
vestigate (3 = 1,2,4,00 to check the dependence on the 
abruptness of the cutoff. 

The entrainment fraction for the winds depends on the 
dumpiness of the IGM, and on instabilities in the outer 
shell of an expanding wind. There is currently no good 
basis - theoretical or observational - for a particular as- 
sumption of ecnt, so we will try values between 1% and 
100% with a fiducial value of 10%. 

The critical SFR/area to drive a wind can be roughly 
addressed using observations. Observational samples 
of 'normal' star-forming disks versus starbursts indicate 
that normal disks tend to have SFR/(area) of 0.001 — 
0.1 M0 yr~^ kpc~ , whereas starburst regions typically 
have 0.1 - lOOOMgyr^i kpc"^ (Kennicutt 1998). Heck- 
man (2000) finds a similar threshold for starburst-driven 
superwinds of 0.1 yr""'^ kpc~ We adopt this value 
as our fiducial SFRcrit- This choice, as will be discussed 
in § leads to winds in most (but not all) galaxies at 
z > 2, and about half of the galaxies at z — 1. However, 
as noted above, the simulated galaxies do not attain the 
very high areal SFRs seen in starburst galaxies because 

shell velocity. 

kpc-2. The reason for this large disparity is unclear. 
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Table 1 
Wind Parameters 



Parameter 


Description 


Default value 


Range 


O-T itcrit 


ciiLicai or n,/ 1 ciisK aieai loi Liieiiiiai wiiiq 


u. iivi0 / yi / Kpc 


n nm n 9 

U.UUi-U.Z 


fid 
^out 


outflow velocity at initial radius 7'o 


600 km/s 


300-1000 


fout 


width of outflow velocity distribution 


200 km/s 




X 


'fraction' of supernova energy in wind (c.f. § 2.2) 


1.0 


0.5-2 


font 


fraction of ambient material entrained in wind 


0.1 


0.01-1 




fraction of galaxy mass enclosed within initial radius 


0.1 


0.05-0.2 




fraction of metals distributed non-locally 


1.0 


0.5-1.0 


J dust 


portion of ejected metals in dust 


0.5 


0.0-0.5 


a 


WS'P cx (r//i)" 


3 


1-4 


p 


Wind and initial shell energy attenuated by (SFR/SFRcrit)^'^- 


2 


l,2,4,oo 



we average the SFR over ^ 10^ yr (longer than a typ- 
ical starburst 'event'), and because the simulations can- 
not resolve the scale of a typical starburst nucleus. Star- 
bursts that should occur in simulated galaxies are therefore 
washed out in both time and space when we compute our 
areal SFRs. Thus we are essentially assuming that sim- 
ulated galaxies exceeding the critical SFR would contain 
the same sorts of starburst regions as observed galaxies, if 
only such regions could be resolved. We also investigate 
'quiescent' models in which we set a smaller SFR/(area) 
of 0.001 M0 yr~^ kpc~^ so that essentially all galaxies in 
the simulation drive weak winds. 

We start the shells at a radius enclosing a fraction ^ 
of the galaxy's baryonic mass, and choose ^ so that the 
initial radii obtained are comparable in scale to the ar- 
eas observed in the studies upon which we base the wind 
mass outflow rates and initial velocities. These observed 
regions are typically 10-100 kpc^; taking ^ ~ 0.1 (our fldu- 
cial value) gives rg ~ 1 — 100 kpc at z = for simulation 
galaxies with circular velocities similar to those of the ob- 
served galaxies. We check the importance of f in § |[ 

Finally, we must choose a form for the distribution func- 
tion W{r,h). The flnal distribution of the wind material 
should presumably be 'piled up' where it is stopped by 
the IGM (or gravity) , as in the bow shocks often observed 
at the end of jets. But any more detailed assumption of 
W{r, h) seems difficult to justify. For simplicity and gener- 
ality, we assume W{r,h) oc (r/h)", for 1 < a < 4 (with a 
default of a = 3), which should span a range of reasonable 
cases. 

2.3. Prescription 3: Radiation Pressure Ejection 

For bright galaxies, the outward radiation pressure felt 
by a dust grain in the interstellar medium (ISM) can ex- 
ceed the inward gravitational pull, suggesting the possibil- 
ity that dust grains can be expelled from galactic disks into 
halos or even into the IGM (see Chiao & Wickramasinghe 
1972; Ferrara et al. 1990; Shustov & Vibe 1995; Davies et 
al. 1998; Simonsen & Hannestad 1999). For a spherical 
model galaxy with a radially increasing mass/luminosity 
ratio Af/i, this outward efflux would inevitably occur 
within some critical radius where the radiation pressure 
and gravitational forces on the grain balance. 



In our prescription to model this physical effect, we dis- 
tribute dust near its progenitor galaxy in a way that re- 
flects the equflibrium distribution of dust 'levitating' at the 
force-balance radius. The maximum radius at which dust 
can be in force equilibrium, ft,'^"''', depends on the galaxy 
luminosity (which we compute using spectral synthesis) 
and on the distribution of mass (given directly by the sim- 
ulation). The radial density proflle of the levitating dust 
d epen ds mostly on the dust properties, described below in 



2.5 



We do not calculate the destruction of grains be fore 
they reach this radius, but we do (as described in § 2.5 ) 



calculate thermal sputtering of dust after it is deposited. 

To implement this prescription numerically, we flrst di- 
vide the gas and star particles into bound groups using 
the SKID algorithm as in the wind prescription, and de- 
termine the center of star formation Tc using Eq.H Then 
about this center we determine a radius /i"^"*^* such that 
the gravitational force due to the matter within /i^ust ]-,g^j_ 
ances the outward radiation pressure due to the group's 
luminosity L^^ i.e. 



1 



dV LuKy — 



GMtotjh) 



(8) 



where Mtot(^) is the total group mass within h, L^, is 
the group luminosity, K,y is the dust cross section/unit 
mass, and we assume spherical symmetry. Assuming a 
bolometric luminosity Lboi and a radiation-pressure effi- 
ciency Qpr(r) (defined as the ratio of Hi, to geometrical 
cross section per unit mass, 7ra^/ {Ana^p/S), averaged over 
a blackbody of a temperature T that reasonably approxi- 
mates the spectrum of a galaxy), we have 



L 



bol 



IGirGcap, 



'dust 



Mtot(/l'l"^*) iQpr 

where a is the grain radius, giving 



Lho\/Lc, 



Mtot/M0(/ld-t) 



1.76 



Qpr(r) 

a[p,m] 



Pdust 

gcm~^ 



(9) 



(10) 



To compute £boi, we track the 'effective' age of each 
star particle, then take Lboi(''") (in units of Lq/Mq) for 
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a stellar population of that age from the stellar synthe- 
sis models of Bruzual & Chariot (1993),^*' for an assumed 
IMF. We then sum this luminosity over all star particles 
in the group: 



L 



bol 



E 



(11) 



To compute r for a newly formed star particle, we assume 
that the SFR is constant between simulation outputs, so 
that the luminosity of stars formed between time steps tp 
and tc can be computed as 



1 



(12) 



We then invert Lboi(''cff) to obtain Toff. When a new stel- 
lar mass Amj is added to an existing star particle of age 
Ti{tp) and mass rriiitp), the new effective age Ti(tc) will be 
given by 



Lho\{n{tc)) 



Lbol{Tes)Ami + Lbol(Ti(tp) +tc- tp)mi 



rrii + Am 

(13) 

As pointed out by Davies et al. (1998), the luminos- 
ity calculated using Eq. |ll| must be corrected for extinc- 
tion by dust in the galaxy, because radiation re-emitted in 
the far-infrared (FIR) provides negligible radiation pres- 
sure. Large dust corrections have been deduced for high 
luminosity local spirals (Wang & Heckman 1996) and star- 
bursts (Heckman et al. 1998), and in high redshift star- 
forming galaxies (e.g., Calzetti & Heckman 1999; Granato 
et al. 2000). As there is significant disagreement as to 
what dust correction is appropriate in any given context, 
the topic requires some discussion. In both starburst nu- 
clei and in disk galaxies, the dust correction (i.e. the ratio 
Lir/Luv) appears to be correlated with bolometric lumi- 
nosity. Heckman et al. (1998), however, find that the 
correlation with metallicity is stronger, so the correlation 
with luminosity (or mass) may be partially the metallicity 
correlation combined with a luminosity-metallicity rela- 
tion. For example, Heckman et al. (1998, Figure 2b) give 
a metallicity-extinction relation of 



log(LiR/i 



uvj 



1.45-1- 1.65 log(Z/Z0), (14) 



where Z is the oxygen abundance relative to solar, Luv 
is luminosity at 0.19 /xm, and Lir is the integrated ~ 
40 — 120 /im luminosity. Combining this with an (ellip- 
tical galaxy) luminosity-metallicity relation 

log(Z/Z0) = 0.4 log(Ls/ Lq) - 4.4 (15) 

(from Zaritsky, Kennicutt & Huchra 1994), we obtain 

log(LiR/Luv) = 0.66 log(iB/ Lq) ~ 5.8, (16) 

that reproduces the slope of the disk galaxy extinction- 
luminosity relation 

log(LiR/Luv) - (O.5±G.l)log(LB/L0)~(3.7±G.7) (17) 



of Wang & Heckman (1996) tolerably well, even though 
Eq. 14 was derived for starbursts, Eq. |l^ for ellipticals, and 
Eq. 17 was derived for spiral galaxies. Neither Eq. nor 
Eq. 17 can straightforwardly be applied as a dust correc- 
tion: they relate the FIR flux to the UV (« 0.19 /^m) radia- 
tion rather than to the UV-optical-NIR radiation (denoted 
henceforth by the subscript 'opt') that drives dust. We 
assume a constant relation between UV and UV-optical- 
NIR attenuation, given by C/ = log(L/fl/Lopt).^^ This 
constant is not 'free', however, as there is the observa- 
tional constraint that the total integrated (over the cos- 
mic history and all galaxies) UV radiation as observed in 
the UV/optical/NIR background is comparable to the to- 
tal far infrared emission in the FIR background radiation. 
Madau & Pozzetti (2000) find that the energy in these two 
radiation backgrounds have a ratio of 1/3 ^ Fib./ F^pt ^ 2, 
with a probable value of ~ 1. Because Eq. |lj along with a 
mass-metallicity relation can yield an extinction-mass re- 
lation, and because the simulations should by themselves 
give a mass-metallicity relation, we adopt 



log(LiR/Lopt) = 1.45 + 1.651og(Z/Z0) + C 



/ 



(18) 



as a fiducial dust correction, with C / set by requiring that 
the total emitted FIR/UV ratio (an output of the cal- 
culations) is near unity. This typically leads to values 
of Cf ^ —1. We also allow for the possibility of a dust 
correction depending purely on luminosity (from Wang & 
Heckman 1996): 



log(LiR/Lopt) = 0.51og(Lboi) - 4.7 Cf. 



(19) 



Here, we use the galaxy bolometric luminosities Lboi hi 
place of Wang & Heckman's Luv + ^fir- Finally, we 
also try a constant dust correction log(LiR/Lopt) = Cf, 
in which case C/ ^ from the background density con- 
straint. 

The dust-corrected luminosity can be used with Eq. |l^ 
to find /i'^"''', the metal mass Aw^^'p is then distributed 
among the gas particles within the radius /i'^"^* with the 
mass profile W^^'^{r, /i'^"*'*) as in the wind prescription. 

For Vl^stp = 5{r — /i'^"''*), the prescription accurately 
captures the physics of dust of a single grain size being 
ejected from a spherical region with sufficient gas drag on 
the dust that the grains do not attain escape velocity. In 
general, the prescription captures the essential physics of 
dust ejection fairly well, but nevertheless has a number of 
important shortcomings: 

1. Charged grains can be confined to galaxies by mag- 
netic fields. Dust might escape if magnetic fields 
have a component perpendicular to the galaxy; 
winds or Parker instabilities exacerbated by the ra- 
diation pressure may enhance such a component 
(Chiao & Wickramasinghe 1972; Ferrara et al. 1991; 
Shustov & Vibe 1995). Ferrara et al. (1991) and 
Davies et al. (1998) have also noted that grains are 
charged only sporadically. Magnetic fields may also 
be much weaker at high redshifts if they have been 
amplified by a dynamo. Determining the effects of 



^''The models are available via anonymous FTP from ftp.noao.edu. 

^^Realistically Cf should be constant only for a fixed effective galaxy temperature; see § ^. Also, since Eq. |l^ applies to starburst regions, 
Cf also absorbs geometrical differences in extinction between starburst regions and the dust corrected region. 
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4. 



magnetic fields would require information about the 
magnetic structure of galaxies, their halos, and the 
IGM that is presently unavailable, so we will assume 
in this study that dust escapes magnetic confinement 
(though we emphasize that it also possible that it 
does not). 

After its deposition, dust is coupled to gas (i.e. fol- 
lows the gas particles). Thus dust distributed in the 
halo can simply fall back into the galaxy, rather than 
maintaining its distance from the galaxy. 

More detailed calculations show that grains often at- 
tain a high velocity at small radii, carrying them past 
^dust g^^^ perhaps to 'infinity'. Thus our estimates 
of the radius to which grains could escape is proba- 
bly an underestimate, because the grains would in- 
evitably reach that radius with some outward veloc- 
ity that would carry them past it. However, if the 
gas drag is high (or if magnetic fields are important), 
the time to escape the galactic disk could exceed the 
dust destruction time in that environment. Thus 
only in certain situations would the dust escape the 
inner galactic region and 'levitate' in the halo as the 
prescription describes. See Aguirre et al. (2001b) 
for some discussion. 

As in the wind case, the method is only fully con- 
sistent if the dust ejection time (i.e. the transport 
time from the disk into its final equilibrium posi- 
tion) is shorter than the interval between simulation 
outputs or between the ejection time and the time 
at which results are desired, and shorter than the 
timescale over which the galaxy properties change. 
The err ors a t late times are unlikely to be important 
(see §§ 2.2, 1^) unless the transport time substan- 
tially exceeds the time between ejection and 'obser- 
vation'. Dust moving with the high velocities (up 
to - lOOOkms-i) found by Shustov & Vibe (1995) 
would not encounter this problem but slower-moving 
dust might. 

The distribution of metals occurs without regard to 
the existence of nearby galaxies whose luminosity 
and gravitation would affect the force balance. Also, 
as in the wind prescription galaxies can (unrealisti- 
cally) enrich each other with their metallic outflows, 
though this is a small effect with the distribution 
method we employ. 



6. 



7. 



8. The dust correction we apply is rather uncertain. 
We have chosen a dust correction to the UV-optical- 
NIR luminosity based on metallicity, and also tried 
a luminosity-dependent correction. The metallic- 
ity correction is reasonable and in accord with ob- 
servations, but it neglects the distribution of dust 
within galaxies, the gas fraction, and also the effec- 
tive temperature of the radiation. At higher redshift, 
galaxies are probably more compact, have higher gas 
fractions, and emit more light at short wavelengths 
where dust attenuation is more effective. All of these 
effects may increase the importance of dust at high 
z more than the decrease in metallicity suppresses 
extinction. But the luminosity-dependent correction 
we try is stronger at high z and can be used as a 
check on our assumptions. 

While our prescription for modeling the ejection of dust 
could be significantly improved given a better understand- 
ing of the dust ejection process and of the correct dust 
correction, it should give a reasonable estimate of which 
galaxies could eject dust, and of the ejection radius. 

2.3.1. Fiducial Parameter Values for Dust Ejection 

The key physical parameters in the dust ejection model 
are the dust absorption efficiency, the dust distribution 
function VF^rp^ t]-^e adjustment factor C/ in the dust cor- 
rection to the galaxy luminosity, and the IMF used in the 
spectral synthesis. 

For the dust absorption efficiency, we adopt the val- 
ues calculated by Laor & Draine (1993) for silicate and 
graphite dust, averaged over a Planck spectrum of T = 
12000 K.^^ We also assume a specific gravity of 2.2 gcm^'^ 
for graphite and 3.3 gcm""^ for silicates. These assump- 
tions yield maximum values of Qpr/ap 19(4) for graphite 
(sificate) grains. For T « 10000 K and T « 8000 K these 
values are 16(3) and 13(2), respectively. These values de- 
termine the maximum radius to which grains could levi- 
tate; but grains of different sizes will levitate to different 
radii, owing to their range in Qpr- 

For a general grain mass-size distribution in mass 
dm{a)/da and absorption efficiency (5pi (a)/a, we can de- 
rive an approximate form of W^'^^{r, h'^^^^)dr by assuming 
a fiat galactic rotation curve, i.e. Mtot(?') oc r. Then the 
force balance equation (Eq. ^ gives r — KQp^{a)/a for 
some proportionality constant if, and we can take 

dm(a) 



W^'^{r) 



da- 



da 



-6{r - KQp,{a)/a) 



(20) 



that should accurately capture the mass distribution of 
levitating grains with a given grain-size distribution. For 
the Qpr values of Laor & Draine with T = 12000 K, and 
the 'FED' g rain size distribution of Kim, Martin & Hendry 
(1994; see § 2.5), the derived W{r, h) can be reasonably fit 
by a second-degree polynomial, which is used in the actual 
algorithm. 

The default value of C/ will be taken as the value neces- 
sary to give Fm/Fopt ~ 1 in the background radiation, and 
the default IMF is Scalo (1986), from 0.1 Mq to 125 M©. 

The fiducial parameter values just described are sum- 
marized in Table 3 and are used in the fiducial models 
listed in Tables |^, |, and |8[ 

^^We do not account for the eflfect of extinction on the effective temperature of the escaping radiation. 



We ignore the radiation-related forces on grains due 
to photodesorption and the photoelectric effect (see 
Weingartner & Draine 2000). Both of these effects 
would increase the radial force on the grain by a 
small factor. Also, we neglect the luminosity due to 
accretion onto black holes. 

Galaxies are not spherical, so the spatial distribu- 
tion of 'levitating' grains will only be approximately 
spherical if h'^^^*^ is much larger than the character- 
istic size of the galaxy. 
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Table 2 

Parameters for Dust Ejection 



Parameter 



Description 



Default value 



Range 



oj 
'dust 

}pr/ap 



correction 
IMF 



fraction of metals distributed non-locally 
portion of ejected metals in dust 

absorption efficiency for dust 
free constant in dust correction 
dust correction type 
IMF for L/M determination 



0.5 

1.0 

19(0), 4(5') 
fixed by im/Lopt 
Z 
Scalo 



0.0-1.0 
0.0-1.0 



Z /Lboi/constant 
Scalo/Salpeter 



2.4. Other Parameters and Considerations 

A few more general considerations (and their associated 
parameters) are common to two or more of our prescrip- 
tions and we discuss them here. The total normalization 
of the metal mass in our calculations is given by the 'effec- 
tive yield' y,, defined as the ratio of metal mass returned 
to the ISM to the stellar mass formed, in solar units. We 
take y*=l (i.e. solar yield) for all runs both because this is 
conventional and because all results can be simply scaled 
to a different effective yield. ^'^ As discussed later, a higher 
overall effective yield may in fact be called for by the obser- 
vations (see also Renzini 1997; Pagel 1999; Aguirre 1999). 

The yield y, is split into gaseous metals and metals 
locked in dust. For metals distributed in the local gas, we 
take the ratio /dust =(dust mass)/(total metal mass)= 0.5, 
as in local galaxies, and as suggested by observations of 
damped Ly-a absorbers (Pei, Fall & Hauser 1999; but see 
Pettini et al. 1997 for a lower estimate). For ejected met- 
als, we allow a different ratio f'^^^^- For wind ejection, we 
should have /^^^^ <, /dust, which is sensible if most of the 
ejected metal mass is in the ISM of the galaxy that has 
been entrained by the wind (it is an upper limit because 
some dust would be destroyed during ejection). The value 
may be somewhat different - in a quite unknown way - 
if most metals are contained in the hot wind itself. For 
radiation pressure ejection, /^^^^ ^ 1 applies; the fraction 
then represents the survival fraction of dust as it traverses 
the halo during its ejection. 

The fraction Foj of metals that are ejected (versus locally 
distributed) could vary anywhere from near zero, for ra- 
diation pressure ejection where gas drag or magnetic con- 
finement is very strong, to near one or more for galactic- 
scale winds where the metal-rich supernova ejecta escapes 
along with some entrained gas. We use a fiducial value 
of = 1 for 'superwinds', and Foj — 0.5 for quiescent 
winds. For dust, Kej > 0.5 could not be maintained for 
much of a galaxy's lifetime since only ^ 1/2 of a typ- 
ical galaxy's metals are in dust at any given time, and 
even values of Yej ~ 0.5 would severely change the abun- 
dance ratios of refractory vs. non-refractory elements in 
the galaxies. Hence, this is probably an upper limit. 

A final parameter is introduced in an effort to cor- 
rect for a possible disparity between the simulations and 
reality: whereas the simulations tend to find Vl^^™{z = 



0) ^ 0.011, values of 0.002 ^ rjf "(z = 0) ^ 0.006 are esti- 
mated from observations (e.g., Fukugita, Hogan & Peebles 
1998). Three possible reasons for this disparity (discussed 
in §Q) are: 1) The simulation fi* is correct but most of the 
'stellar' mass is unobservable, e.g., brown dwarfs, 2) the 
simulations over-estimate the efficiency of star formation, 
or 3) the observations underestimate f2,. In cases 1) and 
2), the simplest reasonable correction to make is to multi- 
ply the SFR and the luminosity of the simulation galaxies 
by = nf%z = 0)/17f"(z = 0) ^ 0.011 w 0.36 wherever 
they are used. The stellar yield y* should also be multi- 
plied by e^,, although if case 2) holds the resulting stellar 
metallicities will be too low. 



2.5. Treatment of Dust and Dust Destruction 

Our method incorporates a fairly detailed treatment of 
IG dust. We track both the dust mass and grain-size dis- 
tribution for each gas particle, and treat the conversion 
of dust to gaseous metals by thermal sputtering. Dust 
is added to pristine gas with a set grain-size distribution, 
but when adding dust to a gas particle that has dust, we 
average the grain size distributions. 

The extinction properties of dust in galaxies appear to 
be well fit by models employing spherical grains with a 
power-law grain-size distribution. The grain sizes neces- 
sary to account for the extinction data range from a,„in ~ 
0.001 fim up to some cutoff (generally either sharp or ex- 
ponential) above ~ 0.2 /im. We represent the grain-size 
distribution by a set of power laws, i.e. 



dN{a) 
da 



cx a 



ak-3 



for flfc < a < flfc+i, k = l..iV, (21) 



for some set of a^, with dN{a)/da a continuous function. 
We implement this by tracking the dust mass dk{i) = 
dmf^^^{ak)/da, where i is the particle number. The dif- 
ferential dust mass values at the selected set of Ok then 
determine the shape of the grain size distribution; clearly 
for N oo this allows for an arbitrary grain-size distribu- 
tion; but in practice the method is 'efhcient' enough that 
a small number (iV ^ 9) provides good accuracy (see the 
upper panel of Figure |l| for a demonstration). 

The total dust mass in a particle is determined by a 

^''The metallicity-dependent dust correction in the radiation-pressure prescription does depend on j/, , but the effect is largely counteracted 
by changing Cf to ensure the proper i^pia/i^opt in the derived optical and FIR backgrounds. 



11 



Table 3 
General Parameters 



Parameter Description Default value Range 



V* 


mean stellar yield (solar units) 


1 






factor multiplying SFR, mout, ^boi 


1.0 


0.355-1 


/dust 


portion of locally-distributed metals in dust 


0.5 





piecewise power-law integration over dk, i.e. 




dust _ dk{i) 



k=l 



/ gfc+i 



(22) 



log(dfc4 

-4 



L/d/c)/log(afe+i/afc) 



For ak - 
replaced by 



where ak - 

1 < e « 10~* a term in the sum is 
dk{i)ak log(afe+i/afe), which is good to O(e^). 

To add dust to a particle, we simply add to the dk{i) 
values of the assumed grain size distribution at a^, such 
that the total mass corresponding to the added distribu- 
tion is equal to the desired additional dust mass. This 
in general will not give a set of dk that yield the correct 
total mass (though it will be close), so we 'renormalize' 
the entire distribution, since the total mass scales with a 
constant multiplying all the dk{i)- This effectively yields 
a 'best fit' of the two summed piece-wise power-law dis- 
tributions using a third with the same mass. A sample 
addition is shown in Fig |l| (lower panel) . 

For dust destruction, we find the sputtered radius Os 
using the gas temperature and density from the simu- 
lation, the time interval since the last simulation out- 
put, and the thermal sputtering yields of Jones et al. 
(1994). The grain distribution is then transformed by 
dN{a)/da dN{a+as)/da, or in practice dk dj,, where 
the latter is given by a power-law interpolation: 



di 



ak 



ak 



ak 



where j is selected so that aj < ak + as < Oj+i . For j > N 
we use j = N — 1, i.e. extrapolation. The fraction of dust 
destroyed could then be calculated by integrating this new 
set of d'f. and comparing to the old. This gives a good ap- 
proximation to the effect of thermal sputtering on the dust 
mass and grain-size distribution. We have not calculated 
non-thermal sputtering of the grains (because it requires 
grain velocities), but it could be important: the grain ve- 
locity Wdust would, in the grain's rest frame, correspond 
roughly to a temperature of 6 x 10^(wdust/100 kms^^)^ K. 



0.001 I 

0.001 




Exoct 

- — — approx 
9-pt, opprox, toto 






' xx - 
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0.001 0.010 0.100 1.000 

Fig. 1. — Demonstration of the dust treatment. Top: Unsput- 
tered (stars) and sputtered (diamonds) grain-size distributions, for 
several values of the sputtered radius, a^, for graphite grains. Dot- 
ted lines show the exact distribution; solid lines show the simulations 
representation using 9 points. Bottom: sputtered distribution with 
as = 0.05 /im added to an unsputtered distribution with half the 
mass, yielding a new grain-size distribution (solid). Exact solutions 
are dotted, 9-point approximations are dashed and solid. 

In this study, for the grain-size distribution we use the 
'power-law exponential decay' (PED) fits of Kim, Mar- 
tin & Hendry (1994) of N(a) oc a" exp(-a/ac). They 
give a — —3.48, ac = 0.28 /.tm for graphite, and a = 
—3.06, flc — 0.14 /xm for silicates. 

While the dust treatment employed here is fairly accu- 
rate, it currently has a few limitations: 

1. The dust destruction is underestimated in the simu- 
lations because we have neglected non-thermal sput- 
tering, and dust destruction during ejection. There- 
fore, we only compute dust destruction correctly af- 
ter the grains have come to rest in the IGM. 

2. Dust destruction may also be underestimated if dust 
encounters very dense, hot regions for shorter time 
intervals than the interval between the simulation 
outputs, for example, this could happen for dust 
cycling through cluster cores. 

3. We currently only treat one grain species at a time, 
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Table 4 

Parameters for Dust Treatment 



Parameter Description Default value Range 



'^min 


minimal grain size 


0.005 Mm 






maximal grain size 


2.0 ^m 




dust type 


dust chemical composition 


graphite 


graphite /silicate 


GSD 


grain-size distribution 


PED 





and add dust to gas with a fixed grain-size distribu- 
tion, even though in a realistic ejection scenario dust 
should be somewhat segregated by size due to differ- 
ences in grain absorption efficiencies and gas drag. 

3. OBSERVATIONS CONCERNING COSMIC METALS 

Sections p3[-[7| will discuss the models we have run, and 
the results obtained. First, however, it will be useful to 
briefly review the existing observations concerning cosmic 
metallicity. 

3.1. Cluster and Group Metallicity 

Measurements of elemental abundances using X-ray line 
emission have revealed that the intracluster medium of 
rich galaxy clusters is highly enriched, to between 1/3 and 
1/2 solar metallicity (e.g., Mushotzky et al. 1996). Be- 
cause the intracluster gas mass in a typical rich cluster 
exceeds the stellar mass in cluster galaxies by a factor of 

5 — 10 (e.g., Renzini 1997), this implies that cluster 
galaxies have probably ejected a large fraction (possibly 
up to three quarters) of their metals. It also seems that 
a super-solar yield is necessary to account for all of the 
metals, since dividing the total cluster metal mass by the 
total stellar mass in the cluster gives a yield of y, ^2 — 4 
(Renzini 1997; see also Aguirre 1999). It should be noted 
that measured cluster metallicities are emission-weighted 
and favor the cluster cores; strong radial metallicity gra- 
dients would imply a different mean metallicity for the 
cluster. Cluster abundance gradients are currently some- 
what inconclusive (see Renzini 2000) , but seem to be weak 
except in clusters with strong cooling flows (White 2000; 
Finoguenov, David & Ponman 1999). 

While there is general agreement that in the cores of 
clusters at z < 0.3 the gas typically has roughly con- 
stant Z/Zq « 0.3 — 0.5 for cluster temperatures down to 
~ 3.0 KeV (e.g., Renzini 1999), the observed enrichment of 
cooler clusters and groups is less secure. Davis, Mulchaey 
& Mushotzky (1999) find that Z/Z - 0.3 for intra-group 
gas persists down to temperatures T ~ 1.5 KeV, but for 
cooler groups the metallicity - and perhaps also the baryon 
fraction - declines, perhaps implying that winds can re- 
move metals even from groups. Buote (2000), however, 
finds that two temperature models fit the X-ray data bet- 
ter than the one-temperature models used by Davis et 
al. and others — these fits yield significantly higher abun- 
dances in groups {Z/Zq ^ 0.5 — 1). Well-resolved observa- 
tions using Chandra should provide much firmer data on 
the enrichment properties of groups. 

^*If a significant fraction of cluster stars are intergalactic, the factor 



3.2. Metals in Ly-a Absorbers 

A very useful window into the chemical enrichment of 
the IGM is provided by studies of quasar absorption lines. 
It is now widely accepted that the 'forest' of Lya lines 
found in the spectra of z J; 2 quasars is due to absorption 
by a smoothly fluctuating neutral hydrogen component of 
the IGM (e.g. Cen et al. 1994; Zhang, Anninos & Norman 
1995; Hernquist et al. 1996). Although these absorbers 
were initially expected to be pristine, high resolution spec- 
troscopy has unambiguously identified metal lines (chiefly 
CIV and SilV) associated with N{H I) ^ lO^^-^cm-^ ab- 
sorbers (e.g., Cowie et al. 1995; Songaila & Cowie 1996; 
Cowie & Songaila 1998; Elhson et al. 2000). Applying an 
ionization correction to the abundances derived from the 
line column densities (relative to the absorbers' H I col- 
umn densities), and using a correlation between the Lya 
column density and gas overdensity from numerical simu- 
lations, these lines can give useful information about the 
metallicity of the low-density component of the IGM. All 
of the absorption line studies essentially agree on an in- 
ferred metalhcity of Z/Zq ~ lO'^-^ for 2.5 ^ z < 3.5 
absorbers with N{H I) ^ lO"'^^'^ cm with about an or- 
der of magnitude scatter in the metallicity for different 
absorbers (e.g., Songaila & Cowie 1996; Ranch, Haehnelt 
& Steinmetz 1997; Songaila 1997; Hellsten et al. 1997; 
Dave et al. 1998). However, the metallicity of lower col- 
umn density regions, corresponding to physical overdensi- 
ties p/p ~ 1, which would represent the bulk of baryonic 
matter at z ~ 3, is more uncertain. The metal lines cor- 
responding to the low-column density Ly-a absorbers are 
generally too weak to detect directly and their presence 
can only be derived in a statistical sense, e.g., by analyz- 
ing the median absorption per pixel (Cowie & Songaila 
1998). Recently, Elhson et al. (2000) used a very high 
quality quasar spectrum to show that the CIV enrichment 
must extend to column densities significantly lower than 
iV(i7/) = 10"-^cm-2. Schayeetal. (2000) demonstrated 
that OVI is a more sensitive probe of the metallicity in 
low-density gas than CIV. Using a pixel analysis, they de- 
tected OVI in gas with t{HT) < 1, which corresponds to 
gas densities around the cosmic mean. (See, also, Hellsten 
et al. 1998; Dave et al. 1998.) 

In summary, the typical metallicity of the high column 
density IGM is 10"'^ - 10^^ solar. Very little is known 
about the variation and the scatter of the metallicity as 
a function of density. Although the presence of metals in 
the low-density IGM has been established using statistical 

might be reduced, but galaxies must still lose ~ 1/2 of their metals. 
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techniques, both the overall mean metallicity of the IGM 
and the fraction of the IGM that is enriched have yet to 
be determined. 

3.3. Metals in Galaxies 

To develop a consistent theory of cosmic metallicity, 
one must take into account not only metals in the IGM, 
but also metals distributed in and near the metal-forming 
galaxies. The chemical properties of galaxies constitutes 
a vast subject (see, e.g., Pagel 1997), and here we will 
only summarize some observations that will be of use in 
evaluating the results of our simulations. 

The abundances of heavy elements in galactic gas and 
stars are known to vary with galaxy mass or luminosity 
and galactic radius, and perhaps galaxy type. All three ef- 
fects are evidenced in the very useful compilation by Zarit- 
sky, Kennicutt & Huchra (1994). They note that nearly 
all observed galaxies demonstrate radial abundance gradi- 
ents, making it difficult to assign a particular metallicity 
to a galaxy. However, by choosing a 'characteristic' radius 
(either an isophotal radius or the disk exponential scale 
length), they can compare abundance properties of various 
galaxies at that radius. This reveals a strong metallicity- 
luminosity relation (see Eq. ^ above and Zaritsky et al.. 
Figure 13) for [O/H] in spiral H II regions. A similar but 
shallower M-Z relation exists in [Fe/H] in stars in ellip- 
ticals (Zaritsky et al. 1994; Kobayashi & Arimoto 1999). 
In both cases the characteristic metallicity can range from 
^1/10 Zq in the smallest galaxies to several times solar in 
the largest. These relations can be meaningfully compared 
to the gaseous and stellar metallicities in the simulation 
galaxies at z = 0, although we do not have information on 
the Hubble type of the simulation galaxies. 

At high redshifts, data concerning the chemical prop- 
erties of galaxies can be gleaned from observations of 
Lyman-break galaxies or from studies of damped Lya 
absorbers. Unfortunately, metallicity information about 
Lyman-break galaxies is extremely limited (Pettini et 
al. 2000) and we cannot make a meaningful compar- 
ison with our simulation metallicities at z = 3. The 
damped Lya absorbers constitute the highest column den- 
sity {N{HI) k, 102°cm-2) features of QSO absorption 
spectra. Abundances have been measured for many of 
these systems, giving metallicities of ~ 1/30 — 1/10 solar 
(e.g., Pettini et al. 1999; Prochaska & Wolfe 2000). It is 
unclear, however, exactly what sort of systems the damped 
Lyman absorbers represent; they may arise from a diverse 
population of objects (e.g., dwarf galaxies, outer disks of 
spirals, etc.; see Pettini et al. 1999 and references therein). 
Since we cannot draw a one-to-one relation between these 
objects and an overdensity or average density of a galaxy, 
we will not compare our simulations to these observations 
in the present study. 

4. GALAXY MASSES AND STAR FORMATION RATES IN 
THE SIMULATIONS 

In this study we apply the method described in § |^ 
to three SPH simulations, performed using the method 
described by Katz, Weinberg & Hernquist (1996). The 
simulations themselves are described more specifically in 
Murali et al. (2001). Briefly, all three use parameters 
f]„ = 0.4 = 1 - 17a, = 0.02/1-2, ^ 0.8 and h = 0.65. 

^^The ESO slice project results of Zucca et al. (1997) are very sin 



The first simulation has 144^^ gas and 144^^ dark particles 
in a (50 /i"^ Mpc)'^ box, giving a dark particle mass of 
6.3 X 10^ Mq and a gas particle mass of 8.5 x 10® M©. The 
other simulations are both (11 /i^^ Mpc)'^ in volume, and 
use 2 X 128^ and 2 x 64^ particles, respectively. The 128^ 
simulation has dark and gas particle masses of 9.8 x 10*" Mq 
and 1.3 x 10'' Mq, respectively. The smaller simulation 
boxes stop at z = 3, and most conclusions regarding the 
universe at z = 3 will be drawn from these. The 144^^ 
simulation runs down to z = 0. The 128"^ and 64 simu- 
lations include an ionizing background of from Haardt & 
Madau (1996) (attenuated by a factor of two to adjust for 
our assumed Sit; see Weinberg, Hernquist & Katz 2000), 
whereas the 144^^ simulation does not. 

For our method to yield useful results, it must be applied 
to simulations that can reasonably reproduce the mass, lu- 
minosity, and SFR distributions of observed galaxies. Here 
we briefly discuss the comparison of the simulation mass 
function and SFR to observations and their role in our pre- 
dictions (for more details see Weinberg et al. 1999, 2000; 
Murali et al. 2001; Katz et al., in preparation. Bullock et 
al., in preparation). The luminosity function, which is gen- 
erated in our calculations of radiation-pressure ejection, is 
discussed in § |^. 

The simulations yield a mass function of galaxies di- 
rectly. An 'observational' mass function can be con- 
structed using an observed luminosity function and a pre- 
scription for converting light to mass. Using either a 
constant M/Lb = 7.5 (in solar units) or a M/Lb func- 
tion derived by Salucci & Persic (1999; Persic & Salucci 
1997), the mass function from the 2dF survey (Folkes 
et al. 1999)^^ survey shows quite good agreement - in 
both shape and normalization - with the 144^^ simula- 
tion for galaxy (gas-|-star) masses between 10^*^'^ Mq and 
2011-75 (Katz et al., in preparation). At lower masses, 
galaxies are unresolved by this simulation: as shown by 
Weinberg et al. (1999), only galaxies with ^ 60 SPH 
particles are well resolved; this corresponds to baryonic 
masses of 5.1 x 10^° Mq in the 144^^ simulations and 
7.9 X 10® Mq in the 128"^ runs. At higher masses the sim- 
ulation exhibits a significant excess of galaxies. This may 
result from inaccuracies in the mass-to- light conversion for 
very massive galaxies, or from incompleteness of the sur- 
veys (e.g., due to surface brightness effects; see Impey & 
Bothun 1997), or due to differences in the way masses of 
real and simulated galaxies are estimated, or due to over- 
production of massive galaxies in the simulations. The last 
two uncertainties are exacerbated by the special environ- 
ments of the largest galaxies, most of which are found in 
cluster cores. 

When integrated, the simulation mass function for the 
1443 runs yields a cosmic density in stars of Sl*(z = 
0) = 0.011. This is somewhat higher than the value of 
= 0.004 ±0.002 derived by Fukugita et al. (1998) from 
observed luminosity functions, or from the integrated SFR 
of Steidel et al. (1999), which yields fi, « 0.006. Likewise, 
the simulation SFRs both at low and high redshift tend to 
exceed the observed values (see Weinberg et al. 1999). A 
similar discrepancy was noted based on earlier simulations 
by Katz et al. (1996) and by Pearce et al. (1999). These 
discrepancies may be caused by an observational underes- 
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timate of ri»(z = 0) and of the SFR (due, for example, 
to surface brightness effects or high production of low- 
mass, dim stars), and/or by simulation galaxies being too 
massive by some roughly constant factor, and/or by the 
overproduction of certain galaxies in the simulations.^^ 

In summary, all of the simulations fail to resolve the 
fraction of the mass function below the simulation reso- 
lution limit. At the same time, the simulations may nev- 
ertheless ouerpredict the number and/or mass of galaxies 
that are resolved. We will address these issues in detail in 
future papers (Katz et al., in preparation; Bullock et al., 
in preparation; see Weinberg et al. 1999 for preliminary 
discussion), but note for the present that the large SFR 
and possible overabundance of galaxies we find could result 
from the single-phase description of the gas in the simu- 
lations, or from the lack the of strong feedback. Feedback 
could suppress star formation both in the wind-driving 
galaxy and in nearby galaxies (see Scannepieco & Broad- 
hurst 2000 for some discussion of the latter effect), but 
there is no good way to model such feedback using our 
current methodology. If enrichment were caused mainly by 
dust ejection, dynamics or metal-rich, quiescent winds, our 
neglect of the winds' effect on the galaxies would be more 
self-consistent, though effects other than strong feedback 
would then be required to suppress any excess of small 
galaxies. 

Any real disagreements between the simulations and re- 
ality will clearly affect the accuracy of the IGM enrich- 
ment calculations of this study. We test a possible way of 
compensating for such disagreements by including a pa- 
rameter that can multiply the SFR (which effectively 
sets the supernova rate in our method), the yield y,, and 
the luminosity of the groups. With = ri°''V^*"S t^^is 
would be an accurate adjustment if the discrepancy were 
due to hidden low mass stars, since most of the radiation, 
metal and energy comes from the most massive stars. If 
simulation galaxies are simply more massive by a constant 
factor 0"™/ri°'^'', the same adjustment would properly ac- 
count for the different galaxy luminosities, galaxy SFRs, 
and overall enrichment of the universe. However, since the 
galaxies would still have their original unadjusted masses, 
we would effectively overestimate the effect of gravitation 
as well as underestimate the SFR per unit mass in the 
galaxies. Finally, if the 'observed' values are underesti- 
mated (and the simulations correct), no adjustment is nec- 
essary; the galaxies that the simulations resolve would be 
like their counterparts in reality. However, as noted above 
it seems likely that the low-mass galaxies (at least) are sig- 
nificantly overproduced in the high resolution simulations 
at z > 3. 

5. RESULTS FOR DYNAMICAL REMOVAL 

The removal of metal-enriched gas from galaxies by dy- 
namical processes has been studied in the context of clus- 
ters (e.g., Fukumoto & Ikeuchi 1996; Abadi, Moore & 
Bower 1999; Quilis, Moore & Bower 2000), and in the cos- 
mological context (Gnedin & Ostriker 1997; Gnedin 1998; 
Gen & Ostriker 1999). Gas may be removed during a close 
interaction between two or more galaxies in which a frac- 
tion of the gas attains escape velocity (see, e.g., Barnes & 

^"The high-resolution runs seem to over-produce small galaxies at z ■ 
(1999). 



Hernquist 1992). Alternatively, loosely bound gas may be 
stripped from a galaxy by the ram pressure of the IGM 
through which the galaxy is moving. Removal of metals 
via mergers or tidal disruption may occur in clusters or in 
the general IGM, whereas ram-pressure stripping is prob- 
ably only important in massive clusters^ 
We have applied the method of 
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to our simula- 
tions to examine the dynamical removal of metals. Metals 
are deposited only in bound groups, thus by calculating 
the IGM enrichment at later times we can directly assess 
the amount of metals that can be removed by dynamics 
alone. The results of these trials are shown in Table || 
and Figs. || and ||. The table gives the overall fraction 
of metals that are outside of bound groups, the metallic- 
ity of the IGM at its mean density and at an overdensity 
S = Pgns/iPgas) = 100, and the mean stellar metallicity in 
galaxies. Unless otherwise noted, we include both gaseous 
metals and dust in computing the metallicities (for dynam- 
ical and wind enrichment, this is at very worst a factor of 
two larger than if dust is not included; for dust ejection 
we show dust and gaseous metal fractions independently). 

Figure || shows the mean metallicity (Zg) of the IGM 
as a function of S, and the mean metallicity of hot (> 
5 X 10^ K) bound cluster/group gas as a function of gas 
temperature. Since the clusters tend to have radially de- 
clining temperature gradients, the plotted cluster temper- 
atures are lower than the core temperatures; we have also 
plotted the mean temperature of the IGM within 100 kpc 
for the hottest five clusters. Note also that the cluster gas 
metallicity is slightly overestimated because some hot gas 
associated with galaxies is included, but an examination of 
radial metallicity profiles indicates that the effect is small. 

The mean metallicity vs. overdensity should be inter- 
preted with care; for a highly inhomogeneous metal dis- 
tribution, the mean can be dominated by a few highly en- 
riched particles and should not necessarily be compared to 
the 'typical' metallicity found in Lya absorbers. The me- 
dian particle metallicity gives a better estimate of typical 
particle metallicities, but since each observed line would 
correspond to gas represented by a number of particles, 
the median metallicity of absorption systems should prob- 
ably be higher than the median for particles (for a highly 
inhomogeneous distribution). Thus we expect that the 
mean and median particle metallicities should lie above 
and below the metallicity that should be compared to 'typ- 
ical' observed metallicities at a given column density (see 
Aguirre et al. 2001b for more discussion). In this pa- 
per we display mainly mean metallicities, which should be 
considered upper limits. 

The results show that dynamical removal is rather inef- 
ficient, removing only a few percent of a typical galaxy's 
metals over its lifetime. The enrichment, as clearly shown 
in Fig H (dot-dashed line), is insufficient to account for 
the metallicity of Ly-a absorbers at z ~ 3, which have 
Z « 10^2.5 down to at least N{HI) « lO"-^ cm'^. 

These results are in substantial agreement with simula- 
tions of individual galaxy interactions, which tend to show 
that only a few percent of an interacting galaxy's mass can 
attain escape energy during an encounter (Barnes 1988; 
Barnes & Hernquist 1992; Hernquist 1992, 1993). Our 

: 3 as compared to the luminosity function computed by Steidel et al. 
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Fig. 2. — Enrichment of the IGM for dynamical removal of met- 
als. Top: Mean metallicity {Z) vs. overdensity for the three simula- 
tions at 2: = 3, and for the 144^ simulation at 2 = 0. Also plotted is 
the 2 = 'D144mt' model as a time-resolution test. The top axis is 
log N{H I) (cm~-^) of an absorber corresponding to the bottom-axis 
overdensity at z = 3, using the Dave et al. (1999a) relation. The 
vertically striped box outlines the approximate current constraints 
from the hya forest at 2 fs: 3 while the horizontally striped box 
shows an extension of these constraints to lower densities (see dis- 
cussion in § W). Bottom: Group/cluster mean metallicity (in hot 
gas) vs. mean temperature of hot bound gas at 2 = 0. Arrows 
show the mean temperature of hot gas in the central 100 kpc for the 
hottest clusters. The vertically striped box outlines the approximate 
constraints from cluster X-ray studies. 

results are somewhat in disagreement with those of Gnedin 
(1998), who found - using high resolution SLH-P'^M simu- 
lations at z ^ 4 - that dynamical enrichment was more 
effective than supernova-driven ejection, and also suffi- 
cient to enrich the low-density IGM. While our results 
support the general idea that some metals escape into the 
z = 3,6 < 100 IGM through dynamics alone, we find 
that the amount is negligible compared to the other mech- 
anisms we consider. We cannot rule out the possibility 
that small (unresolved by our simulations), quickly merg- 
ing galaxies at z ^ 5 could enrich the IGM more than our 
calculations suggest; this would, however, require more ef- 
ficient ejection (or a very super-solar yield at high z), since 
if only a few percent of metals are ejected and mixed thor- 
oughly, at least 10% of gas must form stars with solar yield 
to give a uniform enrichment of ^ 10"^'^ Zq. 

The metallicity of the gas in rich clusters in the 144^^ 
run is an order of magnitude below the observed value 
of « 0.3 - 0.5 Zq (see § |), from which we can conclude 

^^Since our simulations do not provide types for the galaxies, we treat 
as ellipticals when comparing to the stellar metallicities. 



that dynamical removal of metals from massive galaxies 
( 10^" '^ Mq) cannot account for the metallicity observed 
in cluster gas. Since ram-pressure stripping should be most 
efficient in small galaxies, it is important to note that 
wc cannot directly address the importance of pressure- 
stripping of dwarf galaxies in IGM enrichment. However, 
if this process is to account for the observed IGM enrich- 
ment, it must happen at rather high redshift, since the 
mass contribution of M < 10^°-^ Mq galaxies at present 
is quite small, i.e. the enrichment would have to happen at 
a high enough redshift that the galaxy mass function had 
a significant fraction of mass in the small galaxies. Note 
also that we have neglected enrichment from Type la su- 
pernovae in intergalactic cluster stars. We do not expect 
this to be significant in the field because only a few per- 
cent of the simulation stars are intergalactic (having been 
moved into the IGM dynamically), but it may be more im- 
portant in clusters, which may have more IG stars. This 
effect would be best treated by introducing 'delayed en- 
richment' into our simulations. Although inefficient (for 
the galaxies resolved) , dynamical removal is more effective 
in clusters than in the field. Assuming that the ratio of 
cluster stars to cluster gas is comparable to the field value 
of ri*/f2f, ~ 0.23, the ejection fraction /oj.ci in rich clusters 
is 

in the 'D144' model, as compared to 0.04 for all galax- 
ies. (Though if stars form more efficiently in clusters, this 
fraction would be smaller.) 

Panels 1 and 4 of Figure || give the average metallicity 
of stars and gas, respectively, as a function of galaxy mass. 
Interestingly, stellar and gaseous M-Z relations exist even 
though galaxies retain nearly all of their metals. The rela- 
tions occur both because the smallest galaxies have higher 
gas fractions, and because they tend to be younger. These 
can be compared to the plotted lines that are rough fits 
to the mass-metallicity relations found by Zaritsky et al. 
(1994) for ellipticals and spirals spanning ^ 15 B magni- 
tudes.^^ Unfortunately the simulations only have enough 
dynamic range to probe the brightest four magnitudes. 
One should also be cautious about the properties of the 
smallest galaxies, near the resolution limit (vertical line). 
Nevertheless, the 'D144' model does seem to exhibit an 
M-Z relation that is too weak in stars, and the stars have 
significantly higher metallicity than observations indicate. 
More effective feedback would prevent stars from forming 
in high-metallicity regions but cannot cure the problem, 
since the metals would then be present in the gas - but 
the gaseous M-Z relation is quite close to the observations. 
Thus it seems that the observed M-Z relation very likely 
indicates metal ejection beyond pure dynamics. We can 
also see from these figures that although all the results are 
shown for a solar yield, the results cannot be scaled by 
changing by a significant amount without clearly vio- 
lating the observed abundances. Thus while it might be 
argued that dynamical removal of metals could pollute the 
IGM more than we predict given a much higher yield (pre- 
sumably due to an IMF biased toward massive stars) , this 
argument would require the higher yield to apply only at 

all galaxies as spirals when comparing to the gas metallicities, and all 
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Fig. 3. — Top three panels: Mean stellar (top three panels) metallicity of galaxies in the fiducial dynamical, wind, and dust (resp. 
'D144', 'W144' and 'P144') models at z = 0. Bottom axis shows stellar+gas mass, top axis shows S-band magnitude, converted from the 
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high-z, before the bulk of cosmic metals are formed. 

Our results concerning dynamical removal of metals are 
not weakened by the large uncertainties in the assumed 
parameters, since in this prescription only the yield is 
important, and can be constrained within a factor of two. 
However, our predictions also are subject to some numer- 
ical uncertainties. 

First, if the timescales for some dynamical processes 
that remove metals from galaxies were shorter than the 
time interval between the simulation outputs used, our 
prescription might not accurately treat their importance. 
We have checked this by repeating our calculation us- 
ing 27 timesteps (model 'D144mt') rather than 19 (model 
'D144'); the additional steps were chosen to roughly halve 
the interval between snapshots during the epoch when 
most star formation takes place (1 ^ z ^ 3). The dif- 
ferences between these two models are quite small (see 
Figure H), indicating that even with 19 timesteps our cal- 
culations have numerically converged. 

A numerical uncertainty that is harder to address is that 
caused by the limited resolution of the simulations them- 
selves. An accurate assessment of the dynamical removal 
of metallic gas - whether during mergers, through tidal 
disruption, or via ram pressure stripping - depends on the 
ability of the simulations to accurately treat both the IGM 
and the structure of the galaxies. The mean physical inter- 
particle spacing in our simulations is 133(5^^/"^(l-|-z)~^ kpc 
for the 128^ runs and bSAS'^^^^l + z^^^kpc for the 144^ 
runs (where 5 is the gas overdensity), large compared to 
the scale of a typical galaxy. Moreover, the central 200 
kpc of a cluster with 100 ^ (5 ^ 1000 has only 200-2000 
SPH particles; it is therefore doubtful that ram pressure ef- 
fects on cluster galaxies are treated accurately (see Abadi, 
Moore & Bower 1999 for discussion). Unfortunately we 
cannot yet perform resolution tests as we have only one 
simulation complete to z — 0. Our simulations treat 
galaxy-galaxy interactions more accurately but still with 
limitations. Large, thick galaxies are probably represented 
well, whereas low mass galaxies and thin disks will not be 
captured. Therefore, while in principle our type of in- 
vestigation can assess the efficiency of dynamical removal 
quite well, in practice we expect limitations due to limited 
resolution. 

We have attempted to test this effect by comparing the 
efficiency of dynamical metal removal in the 64^ and 128'^ 
runs with the restriction that metals are only added to 
gas particles in galaxies of a fixed mass range. For ex- 
ample, we may add metals only to galaxies with 60 — 120 
(gas-|-star) particles in the 64'^ run and only to galaxies 
with 480 — 960 particles in the 128"^ run, to compare the 
relative efficiency of dynamical removal in galaxies with 
baryonic mass 6 x 10^ — 1.2 x 10^" Mq with different reso- 
lution. In this case, we find that ~ 0.15% of metals are lost 
by well-resolved galaxies in the 64^ run by z = 3, whereas 
« 1.7% are lost in the 128^ run. Curiously, we find that 
as we increase the mass cut for galaxies that receive met- 
als, the ejection fraction increases in the 128'^ run, but 
decreases in the 64'^ run. The difference between the two 
runs can be accounted for either by a resolution effect (i.e. 
galaxies of the same mass lose different amounts of metals 
depending on the number of particles constituting them) , 



or by the difference in the mass function of galaxies (i.e. the 
presence or absence of galaxies small compared to those 
with the metals). Were we to assume that in both runs 
most metals are lost from interactions between galaxies of 
comparable mass, then it would necessarily be a resolu- 
tion effect. But if it were purely a resolution effect, it is 
very difficult to see why the ejection fraction would de- 
crease with the number of particles in the 64^ run while 
increasing in the 128^ run. Thus we suspect that metal- 
loss from the massive galaxies is dominated by interactions 
with lower-mass galaxies in the 128'^ run, but dominated 
by interactions with other well-resolved galaxies in the 64^ 
run (the small galaxies being absent). This dependence on 
the presence or absence of small galaxies makes our reso- 
lution test inconclusive. 



6. MODELS AND RESULTS FOR WINDS 



6.1. Fiducial Wind Model 



As we discussed in § 2.2 there have been numerous in- 



vestigations of galactic winds and their possible role in the 
enrichment of the IGM, and the chief empirical data con- 
cerning this process comes from observations of starburst- 
driven superwinds. Our fiducial wind model, labeled 'W' 
in Table || and in the figures, assumes an outfiow ve- 
locity of 600 ± 200kms~^, with wind efficiency (defined 
by Eq. 0) x = 1, an entrainment fraction ecnt = 0.1, 
and a critical SFR/(area) for driving a strong wind of 
SFRcrit = 0.1 Mq yr^i kpc"^ As described in § these 
are chosen to match 'typical' values derived from the obser- 
vations, where possible. Quantities that are not directly 
observable, such as ecnt, the ejection fraction Y^j, and a 
(which controls the distribution that governs the place- 
ment of metals in the IGM) are given reasonable values 
that we vary in subsequent trials. The fiducial parameter 
values are listed in Table 

We show the results of our fiducial wind model in Ta- 
bles H and 0, and in Figs. |, ^, |^ and ^ We first focus 
on results at z = from the 144^ simulation. As listed 
in the first row of Table ^, the galaxies resolved by the 
144"^ simulation lose about 6% of their metal mass to the 
IGM over their lifetimes. At z = 0, the mean metallic- 
ity of the IGM at its mean density, {Zs=i) = 0.008 Z©, 
{Zs=ioq) = 0.02 Zq, and the mean ICM metallicity for hot 
clusters is 0.04^0. Figure ^ allows the most direct com- 
parison to the principal results from the dynamical model, 
shown earlier in Figure |^, and reveals that winds can en- 
rich the low-density IGM must more effectively than can 
the dynamical removal of metals. 

Figure ^ shows some details as to how this enrichment 
occurs. The first panel plots the wind stalling radius 
''stall as a function of galaxy mass, and shows that large 
galaxies with deep potential wells can retain their met- 
als. Most galaxies resolved by the 144^^ simulation cannot 
drive winds past a few tens of kpc (smaller than the galax- 
ies themselves), though a small fraction can drive metals 
hundreds or thousands of kpc into the IGM, where they 
eventually stall after flowing for up to a few Gyr (see panel 
2). We have plotted all angles for each galaxy, so the ver- 
tical 'stripes' demonstrate the range of radii to which the 
shells propagate in different directions 



'Note also that Pettini et aL (2000) find in their study of the z = 2.73 galaxy MS 1512-cB58 no evidence for a non-standard IMF. 
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Table 5 
Dynamical Models Run 



Model 


Variation 


/iGM 










D144 




0.038 


0.0035 


0.0064 


0.036 


0.87 


D144mt 


more outputs 


0.037 


0.0030 


0.0051 


0.032 


0.87 


D64 


64^ 


0.0025 


1.9 X 10-^t 


4.6 X IQ-^ 




0.63 


D64mt 


64^, more outputs 


0.0023 


3.5 X 10-9t 


4.3 X 10-5 




0.69 


D128 


128^ run 


0.017 


2.4 X 10-5 


0.00035 




0.75 


D12<Slt 


12<S'. k^ss outputs 


0.017 


2.5 X 10-" 


0.00040 




().(il 



Note: all results arc given at 2; = for the 144^ run, and at 2: = 3 for the other two runs, 
fin these runs average metallicities at 5 ;^ 10 should not be trusted, since only a few particles have nonzero metallicity in 
this density range. 
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Fig. 4. Quantities in wind ejection for fiducial 'W model at z = 1.5. Panel 1: Pfiysical shell stalling radius vs. galaxy mass. Panel 2: 
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19 



12 13 14 15 It 



10.0000 
1.0000 
0.1000 

A 

SI 
V 

0.0100 

0.0010 
0.0001 

1.0000 
0.1000 
N 0.0100 

V 

0.0010 



Wind Ejection 




W144, z = 3 

W144, z = 

W12B, z = 3 

W64, z = 3 



2 4 

10g(Pg=s/<P.a,>) 



- Ox 




0.0001 L 
0.5 



1.0 



2.0 

<kT> [KeV] 



W144, z = 0, hot gas o 

4.0 8.0 



Fig. 5. — Enrichment of the IGM for the wind ejection of metals 
in the fiducial 'W model (see Tables hi & H). Plotted quantities are 
as in Fig ^. 

(typically 2-32 angles are used for each galaxy). Panel 
3 histograms the final state of the shells ejected between 
z = 1.75 and z = 1.5. About 2/3 of the shells either stall 
within 2ro or turn around and fall back to within rg. An- 
other third stall after z = 1.5, and a smaller fraction stall 
before z = 1.5 but after reaching 2ro. Finally, a small frac- 
tion are still propagating at z = 0. As shown in the last 
panel, most galaxies at z = 1.5 in the fiducial model are 
assumed to be driving winds (i.e. their SFR exceeds the 
assumed critical SFR), but even galaxies with very high 
areal SFRs may not drive an effective outflow, if they are 
very massive. 

Figure |^ shows the details of the radial integration for 
one angle of one galaxy at z = 1.5. The first panel gives 
the shell, wind, infall and Hubble velocities and the IGM 
sound speed versus the shell radius. The shell, starting 
at ~ 800kms~-^ at rp = 1.4kpc, is quickly decelerated 
by gravity and by the sweeping up of matter, leaves the 
~ 15kpc galaxy at ^ 200 km s""'^ with about 6 times its 
initial mass (see the thin, solid line). Outside of the galaxy, 
the same two factors continue to decelerate the wind (see 
panels 2 and 3), though now the wind ram pressure is also 
important, imparting enough force to keep the shell at 
roughly constant velocity out to several hundred kpc. In 
this example, the wind coasts for a long time, eventually 
stopping after running into a nearby mass concentration 
at ^ 600 kpc after ^ 3 Gyr. 

The winds in the fiducial model enrich the IGM little, 

^^Recall that the mean metallicity is not necessarily comparable to the observed metallicitios if the enrichment is non-uniform; see 
Aguirre et al. {2001b) for discussion. 



bringing the cluster metallicity to ~ 0.04 Zq, not much 
more than dynamics alone (see Table ^ and Figure ||). 
Because the overall metal ejection fraction for winds is 
about 50% higher than for dynamics alone, this indicates 
that the cluster environment is suppressing wind escape. 
More massive groups/clusters are enriched to a rather uni- 
form level, whereas cooler groups show a large scatter in 
their metallicity. This scatter, which may or may not be 
supported by observations (see § ||), contrasts with the 
more uniform enrichment by radiation-pressure ejection 
(as discussed in § |^ below) and persists even in models in 
which the metal ejection efficiency is much higher (as in 
the 'W144max' model described below). 

The first row of Table ^ shows z — i results for the 
fiducial model using the 128'^ run. Here, galaxies have lost 
nearly half of their metals, enriching the S = 1 IGM to a 
mean metallicity of 1% solar, somewhat higher than the 
high end of the observed metallicity of the Lya absorbers 
of similar or greater density (see Fig. ^).^^ Relative to 
z — 0, the greater escape fraction of metals occurs because 
the galaxy mass spectrum is shifted toward smaller mass 
galaxies at high redshift; galactic escape velocities thus be- 
come small compared to the fiducial ('W') model's typical 
outfiow velocity. As discussed in Aguirre et al. (2001b), we 
find that the enrichment of low-density regions is limited 
primarily by the time available for the shells to propagate 
into the IGM. 

The dependence of the metal escape efficiency on galaxy 
mass leads to a slightly steeper stellar M-Z relation than 
does dynamical removal alone, as is evident in Figure ^ 
The second and fourth panels show galaxy stellar and 
gaseous metallicity for the fiducial wind model, and the 
enhanced M-Z relation that wind ejection induces on the 
galaxies. We see a similar but stronger effect in the z = 3 
galaxy metallicities of the 'W128' run (not plotted), but 
there the M-Z relation is evident mainly in the enhanced 
metallicity of ^ 10^° Mq galaxies where the escape speed 
approaches the assumed wind velocity. 

As in the case of dynamical removal, we have checked 
some numerical details of the calculation using a number 
of additional models. First, we have verified that using 
roughly half or twice as many simulation outputs changes 
all results insignificantly (i.e. at a similar level as for the 
same test in the dynamics-only prescription; see Table |^ 
and Fig. ^) . We have checked that all the results are sim- 
ilarly insensitive to the time step used in the integration, 
the radius over which we average when integrating (as long 
as it is not large), and the details of the stalling criterion 
for the shell. The accuracy of the radial integration itself 
has been checked by propagating shells with all forces ex- 
cept gravity turned off; comparing the shell velocity (com- 
puted via integration) to the wind velocity (computed us- 
ing energy conservation as per Eq. |^) tests the integration 
accuracy. 

A few more numerical details make small but noticeable 
differences in our results. Halving or doubling the galaxy 
mass contained within the initial shell radius changes the 
results of our fiducial model only slightly (at worst by a 
factor of two in mean metallicity at the lowest densities). 
The difference between rgtaii and ft-wind, which adjusts for 
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unphysical movement of metals, turns out to have little 
effect; all results of both the 'W128' and the 'W144' mod- 
els are changed by < 10% if this ad just ment is removed. 



A third numerical issue, raised in § 2.2, is that while we 



calculate the SFR by dividing the stellar mass formed in 
a galaxy between time tp and tc by tc — tp, we assume 
that this SFR applies for the entire shell propagation time, 
which may exceed tc — tp. Since the SFRs we compute 
are smooth (i.e. don't depend on short episodes of star- 
formation) they should vary on roughly the Hubble time. 
We have run trials in which the wind is turned off en- 
tirely a Hubble time after its launch, and find unimpor- 
tant changes in the results. This insensitivity would also 
extend to episodic star formation, as long as the episodes 
are spaced more frequently than the timescale for the wind 
to decelerate (so that it does not stall between episodes). 

A final numerical de tail is the number of angles we use, 
Na- As described in 



2.2, we choose Na to ensure ~ 16 



gas particles per angle in the galaxy so that the spac- 
ing between successive particles in radius is smaller than 
the scale over which the physical properties of the shell 
change. This typically results in 2-32 angles per galaxy. 
Using the 144'^ and 128'^ simulations, we have run trials 
using 16, 32, 64 and 128 particles per angle. The results 
are changed very little, indicating that enrichment is quite 
similar whether we use ~ 2 or ~ 16 — 32 angles per galaxy. 
This insensitivity is probably due to the fairly spheroidal 
shape of galaxies and their halos in the simulations, and to 
the fact that winds tend to be either confined, or escape to 



large radii where the distribution radius is limited primar- 
ily by the time constraint. Thus it seems that our results 
are not significantly affected by lack of angular resolution. 

6.2. Other Wind Models 

Having examined the fiducial model, we now turn to 
a set of models in which the simulation parameters have 
been varied. In analyzing these variations we may divide 
their effects on metal distribution into three aspects. First, 
the yield y*, when combined with the simulation's star- 
formation efficiency, determines the total metal content in 
the simulation volume and the rough normalization of the 
AI — Z relation. Second, Wout: SFRcrit, Ccnt, X s-^d de- 
termine the fraction of metals that escape into the IGM; 
Fcj controls this directly, whereas the other three control 
whether outflows occur, whether they are strong enough to 
escape the galaxies and whether the metals get far enough 
away from the galaxy so that they do not fall back. These 
parameters, therefore, directly affect the ratios /igm and 
{Zc\) / {Zgai), as well as affecting the M-Z relation and the 
level of enrichment of the IGM. Third, f^ut, rriout/f^out 
and Eont affect how far metals travel from their progenitor 
galaxies; this is most strongly reflected in (Zg^i) / {Zg^ioq) 
or (Zg^i) / {Zc\) , and in the slope of the curves in Figs. H 
and ^. 

The effects of changing are clear; the output metallici- 
ties of the IGM and galaxies are all proportional to this pa- 
rameter. Comparison of the M-Z relation of well-resolved 
simulation galaxies to the observed relation (see Figure |^) 
indicates ~ 1 is appropriate, and values differing from 
this by more than a factor of '--^ 2 could not account for 
the metallicity of observed galaxies unless /cj ^ 1/2. 

The effects of changing the parameters v^^^_, SFRcrit, 
eent and x are a bit more complicated and we describe 
such variations each in t urn. 

As discussed in § |2.2| , superwinds in nearby starburst 
galaxies appear to be characterized by outflow velocities 
of ^ 600 km s^^, but it may be that faster or slower wind 
velocities better capture the real effect of winds. Models 
'Iv' and 'hv' (see Figure I and Tables I and 0) assume wind 
velocities of 300 ± 100 km s'^ and 1000 ± 200kms-\ re- 
spectively. The high velocity model exhibits the features 
one would expect: more metals escape because the out- 
flows are less easily conflned to small radii, and metals 
travel farther from their progenitor galaxies. Also, the M- 
Z relation becomes somewhat stronger (the mass threshold 
above which galaxies tend to retain their metals increases). 
Both trends are reversed in the low velocity model. Inter- 
estingly, the 'W1441v' model shows that some metals reach 
very low-density regions quite well even if the outflow ve- 
locities only drive metals to < 100 kpc. At z = 3 in the 
high resolution simulations the results are similar. IGM 
enrichment is quite high for the 'W128hv' model, and the 
enrichment would be far higher if not for the time con- 
straint. This is not surprising since the model assumes 
1000 km s^^ winds flowing from dwarf galaxies; it is not 
clear that such winds could be sustained for long with- 
out a catastrophic effect on their hosts. Note, however, 
that the value of SFRcrit is such that the galaxies only 
drive strong winds early on; thus galaxies cannot eject all 
of their metals. The 'W1281v' model is interesting for it 
shows that even low velocity winds can enrich the 6 = 1 
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IGM to 1/1000 Zq (assuming that the overestimation of 
the SFR and fi, by the simulation is not too severe; see 
§ |). We have also tested even lower outflow velocities 
at high z; model 'W12811v' assumes outflow velocities of 
100 ± 50kms~^. Even this model enriches the IGM at 
(5 > 10 to a mean level of ^ 1/1000 Zg. 

As discussed in § p.2[ , galaxies with SFR/ (area) 
< 0.1 Mg yr~^ kpc~^ do not seem to drive observable su- 
perwinds, whereas starburst galaxies that do drive such 
winds can have much higher areal SFRs. But since we 
cannot formulate a rigorous criterion for which galaxies in 
our simulations are driving winds, it is useful to check how 
strongly the enrichment predictions depend on SFRcrit- 
Models 'W1441crit' and 'W144hcrit' (not plotted) reveal 
that halving or doubling SFRcrit changes the results in a 
predictable way. For SFRcrit — 0.05, all galaxies at z = 2 
drive outflows, as do a small fraction of z = galaxies. The 
ejection efficiency is significantly increased, as is {Zs=i)] 
there is a somewhat smaller effect on the IGM metallic- 
ity Further decreasing the critical SFR would have only 
a small effect. These trends are roughly reversed when 
SFRcrit is doubled, although in this case yet higher values 
would lead to progressively more suppression of winds. 

An important uncertainty in our assumptions is the en- 
trainment fraction eont, introduced to account for the fact 
that an expanding wind-driven shell almost certainly frag- 
ments and may either leave a large fraction of itself be- 
hind or significantly reduce its covering factor, or, even if 
fragmentation is not severe, may not sweep up all of the 
ambient medium if the gas is clumped. In principal this 
parameter could be deduced from numerical simulations, 
but current simulations do not follow the shell far enough 
into the IGM or resolve small-scale structure in the gas 
well enough to do this. We have tried values of Ccnt of 1%, 
10% and 100% in model 'le', the fiducial model, and 'he' 
respectively. These reveal that the entrainment fraction 
is quite important (see Figure ^ and Tables ||, |^), espe- 
cially in the lower-mass galaxies (the high mass galaxies 
tend to retain their winds gravitationally) . The differences 
are particularly large in the enrichment of the low-density 
IGM. As in the case of the wind velocity, this sensitivity is 
actually useful, as it could be used to constrain the prop- 
erties of winds given the observed enrichment of the IGM 
(if winds are assumed to be responsible). 

Another parameter that is somewhat uncertain is the 
wind efficiency x that fixes the constant of proportionality 
between the wind energy and the SFR. We have calibrated 
X to reproduce the approximate observed mass outflow 
rates. Models '1%' and 'hx' show the effect of changing this 
efficiency by a factor of two. The effect is significant, show- 
ing that continual driving by the wind is important in the 
shell's propagation. The importance of x decreases with 
increasing Cont (as the shell propagation becomes domi- 
nated by conservation of the initial shell momentum). 

Two final parameters that might be varied (and for 
which a priori values are hard to justify) are a (control- 
ling the steepness of the radial profile of metal distribu- 
tion) and (3 (the sharpness of the energy attenuation be- 
low SFRcrit)- Trials with a = 1,2,3 and 4 have at most 
~ 20% differences in the output quantities listed in Ta- 
ble ^. Different values of 13 give important differences in 
the results, but entirely predictable ones. For /3 = the 



results are similar to the 'Q' model described below (i.e. 
they simulate SFRcrit 0). Choosing larger values of (3 
has no effect for small SFRcrit, a drastic effect if all of 
the galaxies have SFR < SFRcrit, and an effect similar to 
raising SFRcrit itself for intermediate values. 

In summary, we find that the fraction of metals that es- 
cape galaxies, expressed as /igm or (Zci)/(Zgai), is fairly 
sensitive to the assumed outflow velocity and to the en- 
trainment fraction (varies by a factor of 2-3 within our 
assumed range), is slightly less sensitive to the critical 
SFR/ (area) and wind efficiency (varies by less than a fac- 
tor of 2), roughly scales with Y^j, and is insensitive to 
other parameters such as a and f3. Very similarly, the en- 
richment of low-density regions, which requires metals to 
travel farther from their progenitor galaxies, is also sen- 
sitive to Woyt and eont {{Zs=i) varies by a factor of up 
to ^ 10). Low-density enrichment is somewhat less sensi- 
tive to the assumed critical SFR and wind efficiency, again 
roughly scaling with Y^^^ and is not sensitive to a or [3. 
Cosmologically- averaged metal mass is determined by 
but can be constrained strongly by comparing the calcula- 
tions to the observed M-Z relation and cluster metallicities. 

After varying many of the parameters individually, we 
have also varied combinations of parameters to model var- 
ious physically distinct possibilities. First, we have per- 
formed trials with = 0.004/f7"^™(z = 0) = 0.36; this is 
equivalent to dividing SFRcrit, and multiplying x and y*, 
by the same value. As described in § 2.4, e* is introduced 



to account for possible differences between the simulated 
and observed SFR, mass function and fi, of galaxies. Set- 
ting a low accurately mimics an IMF in which most mass 
goes into low-mass stars and brown dwarfs, and somewhat 
less reliably adjusts for over-efficient star formation in the 
simulations. The 'Isfr' models show significantly less ejec- 
tion and a smaller enrichment of the IGM. The smaller 
ejection fraction is largely due to the effectively larger 
SFRcrit- Winds reach smaller radii, due to the lower effec- 
tive value of X- The very low metallicity of the low density 
IGM is a product of the lowered and the lower ejection 
efficiency. The former effect is realistic if il^/rif, is signif- 
icantly overestimated by the simulations. However, if we 
do assume a steep IMF so that many low-mass objects are 
present, we see that > 1 would be required if galaxies 
are to have reasonable metallicities, and this would tend 
to cancel the effect. 



As pointed out in §2.2, we have formulated our model 



to simulate powerful winds from galaxies with the highest 
SFRs. But even galaxies that are not und ergoing violent 
star formation may drive winds (c.f. § 2.2), though these 



must have x ^ 1 most of the time, or they would disas- 
semble the entire galaxy over time (the mass outflow rate 
given by Eq. ^ would always be larger than the SFR). The 
'q' and 'qlv' models were chosen to represent such winds; 
SFRcrit is set low enough so that almost all galaxies drive 
winds. The wind efficiency x is 1/lOth its fiducial value, 
and l^j = 0.5. At low-z and from massive galaxies, IGM 
enrichment in model 'W144q' is almost as effective as in 
model 'W144', both because the wind stopping radius de- 
pends fairly weakly on and the larger fraction of galaxies 
driving winds compensates for the less effective ejection. 
At high-z the 'W128q' model enriched significantly less 
than the fiducial model; this occurs because most galaxies 
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at z > 3 are driving winds even for the fiducial value of 
SFRcrit- The 'qlv' model assumes also that the winds are 
relatively slow-moving. The enrichment in this model is 
very weak at both high and low z, thus if 'quiescent' winds 
are to enrich the IGM significantly, they must have fairly 
high velocity. 

Given the assumption that wind speed does not de- 
pend on galaxy mass, low mass galaxies (mostly unre- 
solved by the 144'^ simulation) should eject metals most 
efficiently. Thus we expect our 144'^ simulations to under- 
predict the IGM enrichment due to winds. Nevertheless 
it is interesting that neither the fiducial model nor any of 
its minor variants can account for the metal enrichment of 
groups/clusters. Because we resolve most of the observed 
z = mass function, this indicates that something other 
than winds provides the bulk of the enrichment, or that 
the enrichment happens at fairly high redshift (where the 
mass function shifts to lower mass galaxies not resolved by 
the 144'^ simulation), or that winds are described by pa- 
rameters somewhat different than those we have assumed. 
To make this point more robust, we have generated a 'max- 
imal' model, 'W144max', which combines a high wind ve- 
locity, high wind efficiency, low entrainment fraction and 
low critical SFR. While galaxies in this model eject ^ 20% 
of their metals, the enrichment of clusters is still only 
about 1/lOth solar, several times smaller than observed. 
Moreover, this is achieved at the cost of a mass-metallicity 
relation significantly steeper than that observed. In fact, 
it does not appear possible to fit both the M-Z relation 
and the cluster metallicity for any set of parameters (i.e. 
without modifying the method). Because of the low res- 
olution of our z = Q simulation it is premature to draw 
strong implications for the enrichment of clusters, but we 
hope to return to this topic in a future study. 

An interesting physical effect we can examine with our 
calculations is the effect of winds on dynamical enrichment: 
if metals are moved into galaxy halos, it seems likely that 
they will be more easily removed by dynamical processes. 
This effect has been seen in detailed simulations by Mu- 
rakami & Babul (1999) of galactic winds in clusters, and 
discussed in the context of 'general' IGM enrichment by 
Ferrara et al. (2000). To investigate this effect we have 
run wind models in the 128'^ and 144^ simulations with 
wind velocities of 300 ±100 km s'^ and 1000 ±200 km s-\ 
respectively. In the first, we generate metals only at z = 5, 
and examine their distribution at 3 < z < 5. In the sec- 
ond, we generate metals only at z = 2, and examine their 
distribution at < z < 2. The results are shown in Fig. ||, 
with the 128"^ runs on the left and the 144'^ on the right. 
The top panels shown median metallicities; mean metal- 
licities are given in the bottom panels. The right panels 
show that at low z the effect exists; some metals deposited 
at (5 > 100 find themselves sX 5 <^ 100 at z = 0. But the 
effect is slight, and is strongest (particularly as evinced by 
the median metallicities) at J ~ 100, suggesting that it 
may be happening primarily in groups and clusters. In- 
deed the left panels show that between z = 5 and 2 = 3 
metals tend to migrate from low- to high-density regions 
(the metals which go 'missing' from the plots are those 
absorbed by stars). Thus it seems that in the 'general' 
IGM, dynamics tend to move metals from moderate den- 
sity regions into galaxies and other higher-density regions. 



rather than distributing them more widely into low-density 
regions. 

6.3. Summary and Discussion 

In summary, we find that many physically distinct 
'types' of winds can enrich the IGM to comparable levels. 
For example, quiescent winds (weak winds from all galax- 
ies) can enrich the IGM at all densities to a similar level 
(within a factor of two) as winds with very strong outffows 
from only the galaxies with the very highest SFRs. The 
predictions are, however, different enough (e.g., in the M-Z 
relation) that with refined simulations and better observa- 
tional data our calculations could be used to constrain the 
physical models of winds themselves. 

While the emphasis of this paper is methodology, we 
can draw some general conclusions regarding wind enrich- 
ment of the IGM using our calculation. We repeat that 
we have not shown that winds can develop in any particu- 
lar type of galaxy (as in e.g., Mac Low & Ferrara [1999] or 
Strickland & Stevens [2000]). Instead, we assume the wind 
properties based on observations, amalgamating the phys- 
ical criteria for a wind into the methodological criterion of 
a high enough areal SFR. Under this assumption, galaxies 
can lose a significant fraction of their metals to the IGM. 
Averaged over all galaxies, we find that typically 5-50% of 
metals are expelled into the IGM, and that the fraction is 
highly degenerate in the assumed parameters. 

Our calculations indicate that galactic winds can enrich 
even the low-density regions of the IGM quite effectively. 
At z = 3, winds escape to large radii, and the ensuing 
enrichment is sufficient to roughly account for the metal- 
licity seen in low-density Lya absorbers (see Aguirre et al. 
2001b for more discussion). Enrichment of the IGM be- 
comes progressively higher, and spreads to progressively 
lower density regions of the IGM, as redshift decreases; 
this is in agreement with the findings of Cen & Ostriker 
(1999). At z = 0, we predict that even quite underdense 
regions of the IGM are enriched to a mean metallicity of 
Z <; 0.005 Zq unless the parameters are pushed to unrea- 
sonable values. 

The main implication we can draw for the enrichment 
of cluster gas is that large galaxies can enrich the gas to 
^ 1/10 of its observed metallicity given our fiducial model 
assumptions about the wind properties. (Some of this en- 
richment is direct, and some probably occurs because met- 
als are moved into galactic halos where they are more eas- 
ily removed by dynamical processes.) This is probably an 
underestimate of the overall importance of winds to cluster 
enrichment because we cannot address the importance of 
low-mass galaxies at high redshift, which would eject met- 
als more efficiently. If winds are to account for the metals 
in clusters, we find that the enrichment must happen at 
relatively high redshift, or that one of our assumptions 
regarding winds must be modified. 

Our simulations also make predictions about the prop- 
erties of galaxies. While our range in galaxy masses is too 
small for a conclusive comparison, we find that the M-Z re- 
lation is nevertheless a good way to break the degeneracy 
between our parameters - especially that between outfiow 
speed and wind efficiency. While there seems to be an 
M-Z relation in observed bright {Mb ^ — 21) galaxies, 
we have difficulty reproducing this relation. Several ef- 
fects may contribute to this difficulty. First, our assump- 
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10g(Pgay<Pgas>) l0g(Pgay<P6as>) 

Fig. 8. — Movement of metals by gas processes after their distribution near galaxies by winds. Only the tinicstcp at z = 5 (left panels) 
or at 2: = 2 (right panels) was enriched. The median (top panels) and mean (bottom panels) metallicity (vs. density) at lower rcdshifts is 
plotted, showing the flows of enriched gas between regions of different density. Metals which seem to 'disappear' with decreasing rcdsliift arc 
those incorporated into stars. We also plot the z = results for local enrichment (no winds) at z = 2 (right panels). This line is missing from 
the top right panel because the median gas metallicity is everywhere zero. 
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tions that all galaxy types have the same outflow veloc- 
ity may be flawed; but experimentation reveals that we 
would have to change this assumption drastically to alter 
the high-mass M-Z relation, partially because the highest 
mass galaxies tend to be in clusters where the ICM sup- 
presses their winds. Second, the brightest galaxies may be 
too extended for our method to treat well (i.e. the winds 
are stopped inside and by the galaxy itself). Third, the 
'average' metallicity determined by observations may in- 
clude fewer of the outlying, low-metallicity stars than in 
the simulations. Finally, it may be that something besides 
wind ejection causes the observed M-Z relation in massive 
galaxies. 

The results reported in this section point to a few gen- 
eral conclusions regarding our methodology. First, like 
semi-analytic theories of galaxy formation, the method 
employs a number of parameters that strongly affect the 
calculation predictions and that cannot be ascribed defi- 
nite values using observational data. This limits the num- 
ber of specific general conclusions we can draw. Never- 
theless, the method is excellent for generating predictions 
given a specific model (e.g., the model based on locally 
observed starburst-driven superwinds). Further studies, 
either observational or theoretical, producing more robust 
physical models of winds will generate correspondingly 
more robust predictions of the cosmic metal distribution 
using our method. 

7. MODELS AND RESULTS FOR DUST EJECTION 
7.1. Fiducial Dust Model 

Radiation-pressure driven efflux of dust from galaxies 
provides an interesting alternative to galactic winds to 
pollute the IGM with metals. A number of studies have 
pointed out that typical spiral galaxies could eject a sig- 
nificant amount of their dust if the dust can decouple from 
bound gas and magnetic fields (e.g., Chiao & Wickramas- 
inghc 1972; Ferrara et al. 1990; Shustov & Vibe 1995; 
Davies et al. 1998; Simonsen & Hannestad 1999). Our 
method does not directly address the question of whether 
dust can decouple but can yield a reasonable estimate of 
the equilibrium radius of dust for galaxies with proper- 
ties given by the simulations. If dust can escape the inner 
galaxy there should be no obstacle to its reaching the equi- 
librium radius, provided it can do so in a short enough 
time. We can, therefore, give plausible estimates of the 
radii to which dust grains might be ejected by galaxies of 
various masses, ages and types, at various redshifts, and 
track the distribution of metals after they are deposited, 
even if they only reach the halos of galaxies. 

Our fiducial dust-ejection model, labeled 'P144' and 
'P128' in Table H and in the figures, assumes graphite 
grains with Qpi/ap — 19 (see eq. |l^); this is the maxi- 
mum absorption efficiency from graphite grains subject to 
a 12,000 K blackbody spectrum. We assume a Scalo IMF 
in calculating the stellar luminosities, and a dust correc- 
tion depending on metallicity, adjusted so that the inte- 
grated background light has equal parts in FIR and UV- 
optical-NIR. We assume that half of each galaxy's metals 
are distributed non-locally (Fej = 0.5), which is an up- 
per limit since only about half of a typical galaxy's metal 
mass can be in dust. (Also, the connection between the 
metal outflow rate and the SFR is much less clear for dust 



ejection than for supernova winds.) 

The results of the flducial model are given in Table ^ and 
Figs. ^, ^, [l^, and |ll|. The table and the flrst two flgures 
give the same quantities that were presented for winds, 
while the last two figures give details about the ejection 
from individual galaxies. 

Under our assumptions, dust ejection is fairly efficient. 
In the 144^ simulation, z = galaxies have lost 16% of 
their metals, meaning that 32% of the dust that we as- 
sume can escape actually does escape, with the rest being 
retained either because it falls back into galaxies or be- 
cause the galaxies have insufficient light /mass ratios to 
eject it. The ejected dust enriches the intracluster gas to 
~ 1/5 Zq , and the z — IGM at mean density is enriched 
to 0.02 Zq. These values can be roughly hnearly scaled to 
lower values of Ycj. Dust ejection at high z pollutes the 
(5 = 1 IGM at z = 3 to ~ 10~^ Zq. The values given are 
for all metals (solid and gaseous) . Figure |l^ also shows the 
amount of dust converted to gaseous metals by z = 3 due 
to thermal sputtering, in a number of models including 
the fiducial model. Nonthermal sputtering and dust de- 
struction during ejection would also destroy some dust, so 
the actual level of observable gaseous metals should lie be- 
tween the 'gaseous metals-|-dust' curves and the 'gaseous 
metals only' curves. 

A few aspects of the results for the fiducial model merit 
attention. First, the enrichment of low- mass groups is 
much more uniform than in the wind or dynamics-only 
model (compare figures ||, || and ||). For dust ejection, 
the metallicity of hot gas declines steadily with gas tem- 
perature, leveling off for the hottest clusters. For wind 
ejection - even for runs with /oj ~ 1/5 - there is a large 
scatter in the metallicity of groups. It is currently un- 
clear which of these cases has more observational support. 
Second, as is evident from Figure |[ the M-Z relation of 
galaxies is enhanced by dust ejection over the dynamical 
removal case and is more similar to the observed relation 
(the gas metallicities match the observations particularly 
well). The M-Z relation is, however, somewhat different 
than for winds; although the general trend is similar, there 
is more scatter of low-mass galaxies into high metallicities, 
and a somewhat less abrupt drop at masses M <^ 10^^ M© 
(especially compared to wind models with a high escape 
fraction). As in the wind case, the stellar M-Z relation for 
the most massive galaxies is flatter than observed. 

We can investigate the enrichment process in more de- 
tail by examining the properties of the galaxies driving the 
dust outflows. Figures |l^ and [ill show the relations be- 
tween galaxy mass M, metallicity Z, mass-to-light ratios 
M/ Lb and M/Lboi, and maximum radius of dust ejection 
^dust Interestingly, at z ^ 1, fairly large galaxies drive 
dust to the largest radii, whereas at z ^ 1, dust ejection 
is most efficient in the smallest and the largest galaxies 
(see panel 4 of both flgures). This occurs because large 
galaxies have a relatively larger ratio of stars to gas and 
dark matter. At lower redshifts, this effect is overwhelmed 
in the smallest galaxies by the dust correction (which in- 
creases M/L more in large galaxies; see panels 2 and 5), 
and by the tendency of smaller galaxies to be younger. 

At high z, graphite dust can be driven to ^ IMpc by 
at least some galaxies, whereas at low z galaxies can drive 
dust only to a few hundred kpc. The z = galaxies have 
typical bolometric (dust corrected) M/L values 
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Table 6 
Wind Models, 144^ simulation 



Model 


Variation 


/iGM 




(.^5=100) 


(^cl) 


(■^gal) 


W 144 


none 


U.UOY 


O.UUo 


o.uty 


A OtOT 
0.U3/ 


A 


W i44ne 


Cent — J- 


U.U44 


U.UUoz 


A A1 


A HQ 1 




W1441e 


Cent = 0.01 


0.086 


0.013 


0.035 


0.054 


0.79 


W1441sfr 


= 0.355 


0.035 


0.0015 


0.0025 


0.0070 


0.31 


W144„p32 


32 galactic gas particles/angle 


0.056 


0.0073 


0.019 


0.036 


0.83 


W144„pi28 


128 galactic gas particles/angle 
u^^t = lOOOkms-i 


0.056 


0.0073 


0.018 


0.036 


0.83 


W144hv 


0.088 


0.015 


0.035 


0.057 


0.79 


W1441v 


^ 300kms-i 


0.038 


0.0040 


0.0091 


0.025 


0.85 


W144hcrit 


SFR„it = 0.05 


0.039 


0.0048 


0.0088 


0.023 


0.86 


W1441crit 


SFRcrit =0.2 


0.091 


0.012 


0.039 


0.051 


0.79 


W144hx 


X = 2.0 


0.073 


0.011 


0.028 


0.046 


0.81 


W1441x 


X = 0.5 


0.047 


0.0057 


0.014 


0.032 


0.84 


W144inax 


hv, Icrit, hx, le 


0.22 


0.042 


0.10 


0.12 


0.66 


W144q 


SFR^Ht = 0.001, X = 0.1 


0.052 


0.0048 


0.017 


0.032 


0.83 


W144qlv 


q+lv 


0.036 


0.0039 


0.0082 


0.022 


0.85 



Note: all results are given at z = Q. 



Table 7 
Wind Models, 128^ simulation 



Model 


Variation 


/iGM 


{Z8=l) 


(.^5=100) 




W128 


none 


0.47 


0.01 


0.18 


0.34 


W128he 


Cent = 1 


0.19 


0.00079 


0.09 


0.52 


W1281e 


Cent = 0.01 


0.58 


0.022 


0.19 


0.26 


W1281sfr 


e, = 0.355 


0.18 


0.0017 


0.021 


0.18 


W128„p32 


32 galactic gas particlcs/anglc 


0.47 


0.0093 


0.19 


0.33 


W128npi28 


128 galactic gas particles/angle 
= lOOOkms-i 


0.54 


0.0099 


0.23 


0.28 


WI28I1V 


0.57 


0.017 


0.20 


0.27 


W1281v 


v^^^ = 300kms-i 


0.22 


0.0017 


0.098 


0.51 


W12811V 


v^^^ = lOOkms-i 


0.086 


0.00016 


0.038 


0.61 


W128hcrit 


SFRcrit - 0.05 


0.32 


0.0080 


0.12 


0.40 


W1281crit 


SFRcrit =0.2 


0.48 


0.010 


0.19 


0.33 


W128hx 


X = 2.0 


0.56 


0.016 


0.21 


0.28 


W128Ix 


X = 0.5 


0.36 


0.0056 


0.15 


0.41 


W128q 


SFRcrit = 0.001, X = 0.1 


0.083 


0.00031 


0.036 


0.63 


W128qlv 


q+lv 


0.033 


5.8 X 10"^ 


0.010 


0.68 



Note: all results are given at = 3. 
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Fig. 9. — Enrichment of tlie IGM for radiatiotL-presspre ejection 
of dust in tlie fiducial 'E144' model (see Tables H and kI). Plotted 
quantities are as in Fig Q. 

of ^ 1 — 5 and uncorrected M/Lb ratios of 1 — 4. Panel 
6 of Figure |ri| shows the B-band luminosity function of 
the z = Q simulation galaxies computed using the simple 
spectral synthesis described in 2.3, and compares it to the 
2dF observed luminosity function (dashed line). The sim- 
ulated luminosity function with no dust correction (light 
solid line) is a poor fit to the observations. The dark solid 
line shows the same data with the metallicity-depcndent 
dust correction (given in Eq. used in calculating dust 
ejection. Applying this correction requires a choice of C/ , 
where 10'-^^ (in this context) is the ratio of B-band extinc- 
tion to UV (1900 A) extinction. The curve shown assumes 
C f — —0.5; with such a choice the simulations can roughly 
fit the observed luminosity function. The actual relation 
between UV and blue extinction depends not only on the 
dust absorption curve but also on the amount and distri- 
bution of dust in the galaxy (see, e.g.. Chariot & Fall 2000) 
and is uncertain; but C/ = —0.5 is not an unreasonable 
value, indicating that the dust correction and spectral syn- 
thesis methods can probably produce satisfactory (good 
to within a factor of a few) luminosities for our simulation 
galaxies. 

7.2. Other Dust Models 

We now investigate the effects of variations in the model 
parameters. As in the wind models, determines the to- 
tal simulation metal mass. The amount of metals driven 
from galaxies depends directly on Y^, less directly on 
Qpr/ap, the M/ L ratio from stellar synthesis, and the dust 
correction, and quite weakly on the form of Wgrp ; the last 



four quantities also determine how far the dust travels. In 
turn, Qpr/ap depends primarily on the dust type, M/L 
values depend on the IMF, and the dust correction de- 
pends on the FIR/optical ratio of the extragalactic back- 
ground light. 

The absorption efficiency of dust grains depends mostly 
on their composition; in this study we have used opacities 
for both graphite grains and less efficiently absorbing sili- 
cate grains. Model 'sil' assumes the latter. Since the maxi- 
mal Qp-c/ap is 1/5 that of graphite, ejection is significantly 
less efficient (see Table ^ and Figure |l2|). This is impor- 
tant because the majority of dust mass is locked in silicates 
in most carbon/silicate grain models (e.g., Weingartner 
& Draine 1999; Duley, Jones & WiUiams 1989; Mathis 
& Whiffen 1989). Thus a model in between the fiducial 
model and the 'sil' model is a more accurate representa- 
tion of a realistic dust distribution for a two-component 
model. 

We next tried a model with an assumed Salpeter IMF, 
with 0.1 Mq < M < 125 Mq. This resulted in somewhat 
brighter galaxies at z ^ 0.5, and somewhat dimmer galax- 
ies at z 0.5; the overall effect is to slightly enhance dust 
ejection. The effect on the M-Z relation is also weak. The 
simulation results are also somewhat sensitive to the choice 
of the low-mass cutoff in the IMF, as this changes the stel- 
lar mass with little effect on the luminosity. For example, 
starting the Salpeter IMF at 0.2 Mq rather than 0.1 Mq 
would lower M/L by « 24% and the low- mass IMF advo- 
cated by Gould, Flynn & Bahcall (1996) would decrease it 
even more. These uncertainties are, however, unlikely to 
change M/L by more than a factor of two and are hence 
contained within the range of variations we try. 

The assumed dust correction is more important than the 
chosen IMF, and comparable in importance to the choice of 
grain properties. Maintaining the dust correction based on 
metallicity, we have varied C/ to reproduce a present-day 
ratio of 1/3 < Frm/^opt < 2 in the z = extragalactic 
background (where i^opt includes UV and NIR light also). 
The extremes of this range are shown in models 'wcorr' 
and 'scorr' in Table ||. Although h'^^^*' is generally sev- 
eral times higher in the 'wcorr' model than in 'scorr', the 
enrichment of the low-density IGM is quite similar. This 
indicates that very low-density IGM regions can be spa- 
tially close to high-density metal-forming regions. We also 
modeled a constant dust correction with i^FiR/-Fopt = 1, 
in model 'ccorr', and a luminosity-dependent dust correc- 
tion (model 'Icorr'). The results are only slightly different 
from the fiducial case. As in the other prescriptions, the 
dust ejection calculations are insensitive to the number of 
time steps employed. 

Ejection of dust by radiation pressure requires that the 
dust decouple from the gas. Even in the brightest galaxies, 
radiation pressure cannot overcome the gravitational force 
on both dust and an associated gas mass more than 100 
times as large. We have verified this in our simulations 
using a trial with Qp-c/ap = 0.19, and no dust correction. 
The results are nearly identical to the dynamical removal 
run, indicating that radiation pressure is in this case inef- 
fective even at moving dust into the halos of galaxies. 

7.3. Summary and Discussion 

The fiducial dust ejection model can account for both 
the low-density IGM metallicity at z = 3 and a signifi- 



28 




Fig. 10. — Quantities used in dust ejection for the fiducial 'P144' model at z = 2. Panel one: Bolometric mass/light ratio M/Lboi vs. 
galaxy mass. Panel two: M/Lbol vs. galaxy average metallicity. Panel three: B-band M/Lg vs. galaxy mass. Panel four: Maximal 
dust ejection radius /jdust .^g galaxy mass. Panel five: ftdust mean metallicity. Panel six: B-band luminosity function of simulated 
galaxies, uncorrected for dust (solid, thin) and corrected for dust using metallicity (solid, thick). Top axis gives B-band magnitude. 
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Fig. 11. — Quantities used in dust ejection for the fiducial 'P144' model at z = 0. Plotted quantities are as in Figure |ll| except that we 
include the observed 2dF B-band luminosity function in panel 6 for comparison with the simulation. 
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at z = 3 (right panels). Plotted quantities are as in the top panel of Fig eI Top panels show total (gas+dust) abundances and bottom panels 
show only gas phase abundances (after thermal sputtering). 
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cant fraction of cluster metals. It also fits the stellar M-Z 
relation reasonably well at intermediate masses, but not 
at high masses (where none of our models predict high 
enough metallicity). 

Dust can be ejected from galaxies by radiation pressure 
if three basic conditions hold. First, the luminosity/mass 
ratio toward the galaxy must be large enough that radia- 
tion pressure exceeds the gravitational force. Second, the 
gas drag on dust must be small enough that dust can pass 
through the dense gaseous part of the galaxy and into the 
halo/IGM before being destroyed. Third, the dust must 
not be confined by magnetic fields that are bound to the 
gas. In this study, we have addressed only the first con- 
dition, distributing the dust at the radius where gravita- 
tional and radiation pressure forces balance. We have not 
addressed whether dust can truly decouple from the gas 
or, if so, what dust outflow rate ensues. 

The trials presented in this section indicate that with 
reasonable assumptions about the grain properties, spec- 
tral synthesis and dust correction and with neglect of ex- 
tra light from quasars and extra force on grains from the 
photoelectric effect - each of which gives approximately a 
factor of two uncertainty in the radiation pressure force - 
a significant fraction of the dust can be expelled from the 
simulated galaxies. This conclusion would change only if 
several of the uncertainties conspired to make the force 
an order of magnitude or more smaller than we have esti- 
mated. 

The main weakness of our method is a poor understand- 
ing of dust ejection itself, and in particular how well mag- 
netic fields can confine grains to galaxies. Assuming that 
grains can escape, our method can provide strong predic- 
tions of the spatial distribution, grain-size distribution and 
destruction of dust, as well as the extinction and redden- 
ing due to IG dust. These predictions will be presented in 
detail in a future study. 

8. CONCLUSIONS 

We have developed a new method of calculating the 
chemical evolution of galaxies and of the IGM, with par- 
ticular emphasis on the physical mechanisms that remove 
metals from galaxies. The method is applied to already- 
completed cosmological simulations, so it can be run 
quickly to test changes in the assumed prescriptions and 
parameters. In the method, metals are instantaneously 
placed in gas that is in and nearby galaxies, according to 
parameterized physical prescriptions that estimate where 
the metals would 'land' in the gas after ^ 10^ — 10^ yr. 
We have discussed in detail the prescriptions used in this 
study, pointing out which physical effects are captured by 
them. The method can be used to predict the cosmolog- 
ical distribution of metals given assumptions about their 
ejection, or conversely to study the ejection assumptions 
by comparing our calculations to observations in detail. 

In this paper we have applied the method to several 
cosmological simulations. Using these results we can draw 
some conclusions both about the methodology we have de- 
veloped, and about metal ejection in the real universe. 

In our simulations, removal of metals by purely dynam- 
ical processes such as ram-pressure stripping or tidal dis- 
ruption of galaxies, is relatively ineffective at polluting the 
IGM. Averaged over their mass function, galaxies of mass 
^ lO^*^-^ Mq (i.e. those resolved by our 144^ particle simu- 



lation) lose ^ 4% of their metals by z = 0. This accounts 
for only about l/12th of the metal density observed in the 
gas of rich clusters, though we cannot address the dynam- 
ical enrichment of clusters by smaller galaxies that could 
have comprised a significant fraction of the mass function 
at high z. Dynamical removal alone also cannot account 
for the mass-metallicity relation observed in present-day 
galaxies of all masses. Galaxies of mass M ^ 10^-^ Mq 
also lose ~ 2% of their metals by 2 = 3, enriching Lya ab- 
sorbers with NiHI) - 10"-^cm2 to - lO"'^-^ Zq, about 
100 times less than observed. These results indicate that if 
dynamical removal were to account for metals in the low- 
density Lya absorbers (Gnedin & Ostriker 1997; Gnedin 
1998), it must have been by smaller galaxies, presumably 
at very high rcdshifts (z ^ 6). We also find that the dy- 
namical removal of metals is enhanced if the metals are 
moved into galaxy halos (i.e., by winds or dust-ejection), 
but only slightly and only at low z; at high z the net effect 
of dynamics is to move metals from low- to high-density 
regions. 

If the metals in the IGM came from fairly massive galax- 
ies at z ^ 6, some mechanism other than dynamical re- 
moval must have played an important role. Supernova- 
driven winds are a plausible candidate. Our prescription 
assumes that winds develop at a critical SFR/area, with a 
fixed velocity, and an energy in the wind proportional to 
the SFR. We find that the degree of IGM enrichment is 
not very sensitive to the wind efficiency (unless it is very 
different than we have assumed) or how the metals are 
distributed within the wind 'stopping radius'. The results 
are sensitive to the assumed fraction of the ambient mate- 
rial entrained by the wind, the wind outflow velocity, and 
the critical SFR. If the latter is chosen so that a signifi- 
cant fraction of high-z galaxies drive winds (as indicated 
by observations of Lyman-break galaxies; see Steidel et al. 
2000 and Pettini et al. 2001), then winds with outfiow ve- 
locities of ~ 200kms~^ or more (as also indicated by the 
observations) can escape; to large; <:lista.iic:c;s. and enrich the 
low-density IGM to roughly the level observed at z ~ 3. 
Whether the enrichment can match the observed metallic- 
ities in detail, and whether the wind process itself would 
disturb the low-density IGM more than allowed by obser- 
vations is an important and open question (see Theuns, 
Mo & Schaye [2001] for some discussion of the latter). 

At lower redshifts, winds even from massive galaxies 
may be important, though in our models metal ejection 
from galaxies of M ^ 10^°'^ Mq probably cannot account 
for all of the metals in cluster gas unless a rather extreme 
model is adopted or one of our methodological assump- 
tions is changed. Since smaller galaxies eject metals more 
efficiently, it is possible that winds could account for the 
ICM metallicity if these were included. Wind ejection also 
leads to a mass-metallicity relation comparable to - but 
somewhat steeper than - that observed. Using future sim- 
ulations with a larger range of galaxy masses, the M-Z re- 
lation should be a useful diagnostic of outflows, as it helps 
break the degeneracy between our model parameters in 
determining the metal ejection efiiciency. 

The ejection of dust by radiation pressure is another in- 
teresting way that metals may escape galaxies. We assume 
that a significant fraction of galactic dust escapes to the 
radius where gravitational and radiation-pressure forces 
balance. We find that our basic results are not strongly 
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Table 8 

Radiation Pressure Models Run 



iVioaei 


• • 

Variation 


/IGM 


(^5=l) 


\^5=100/ 






<--/ 


P144 


none 


0.16 


0.018 


0.087 


0.19 


0.71 


-0.96 


P64 


64'^ run 


0.32 


0.00072 


0.08 




0.39 


-0.96 


P128 


128^ run 


0.27 


0.00089 


0.14 




0.48 


-0.96 


P144salp 


Salpeter IMF 


0.19 


0.019 


0.11 


0.20 


0.68 


-0.81 


P128salp 


SALP, 128^ 


0.33 


0.0028 


0.15 




0.45 


-0.81 


P144sil 




0.04 


0.005 


0.0084 


0.03 


0.86 


-1.05 


PI 28sil 


STT; 128-^ 


0.033 


n 00028 


n 0055 




0.65 


-1.05 


P1441sfr 


e* = 0.355 


0.061 


0.003 


0.0065 


0.021 


0.30 


-0.33 


P1281sfr 


= 0.355, 1283 


0.098 


0.00018 


0.012 




0.20 


-0.33 


P144scorr 


-Fprn/^opt = 2 


0.11 


0.015 


0.052 


0.12 


0.77 


-0.64 


P128scorr 


sCORR, 128^ 


0.24 


0.00064 


0.12 




0.50 


-0.64 


P144wcorr 


^FrnZ-Fopt = 1/3 


0.20 


0.019 


0.11 


0.22 


0.67 


-1.26 


P128wcorr 


^FrnZ-Fopt = 1/3 


0.29 


0.0012 


0.15 




0.47 


-1.26 


P1441corr 


Luminosity correction 


0.15 


0.016 


0.073 


0.16 


0.73 


-1.09 


P1281corr 


Luminosity correction 


0.28 


0.0011 


0.14 




0.48 


-1.09 


P144ccorr 


const, correction 


0.15 


0.016 


0.077 


0.18 


0.73 





P128ccorr 


const, correction 


0.22 


0.00054 


0.12 




0.51 







Note: all results 


are given at z = for the 144^^ run, 


and at z = 


3 for the other two 


runs. 







sensitive to the assumed IMF or dust correction. The re- 
sults do, however, depend on the grain type and on the 
very important assumption that the grains are not con- 
fined by the gas or magnetic fields in galaxies. The results 
show that metals removed from ^ 10^'^ Mq galaxies as 
dust, then destroyed in the IGM or galaxy halos by sput- 
tering, could also account for the mean level of IGM en- 
richment observed at z = 3 - although again it is unclear 
whether the distribution agrees with the observations in 
detail. At low redshift, massive galaxies can enrich the 
IGM to the observed levels. Moreover, since dust ejection 
does not 'avoid' high-pressure regions as winds do, dust 
can enrich intragroup gas more uniformly than winds. En- 
richment of the IGM by dust would provide a number of 
chemical signatures; in particular, non-depleted elements 
such at N and Zn should be underrepresented in the IGM 
(see Aguirre et al. 2001b). 

Most generally, our simulations support the view that a 
significant fraction of cosmic metals lie in the IGM, and 
our method provides a useful way to generate predictions 



of the cosmic distribution of metals usable in a number 
of ways. With higher-resolution simulations we should be 
able to more effectively test the importance of low mass 
galaxies at all redshifts, as well as perform more careful 
resolution tests of our results. Future observations, as well 
as more detailed small-scale simulations, will help to de- 
velop more accurate ejection prescriptions. 

This work was supported by NASA Astrophysical The- 
ory Grants NAG5-3922, NAG5-3820, and NAG5-3111, by 
NASA Long-Term Space Astrophysics Grant NAG5-3525, 
and by the NSF under grants ASC93-18185, AGI96-19019, 
and AST-9802568. JG was supported by NASA Grant 
NGT5-50078 for the duration of this work, and AA was 
supported in part by the National Science Foundation 
grant no. PHY-9507695 and by a grant in aid from the 
W.M. Keck Foundation. The simulations were performed 
at the San Diego Supercomputer Genter, and the post- 
simulation processing used a computer purchased under 
the National Science Foundation grant no. PHY-9507695. 



REFERENCES 



Abadi, M. G., Moore, B. & Bower, R. G. 1999, MNRAS, 308, 947 
Abel, T., Anninos, P., Norman, M. L. & Zhang, Y. 1998, ApJ, 508, 
518 

Aguirre, A., Hernquist,.Ti , Weinberg, D „ Katz, N., & Gardner, J. 



2001a. ApJ, in press; [astro- ph/0006345 
Aguirre, A., Hernquist, 



dner, J., & Weinberg, D. 



2001b. ApJ, in press; istro-ph/0006346 

Aguirre, A. 1999, ApJ, 52b, 583 

Balsara, D., Livio, M., & O'Dea, C. P. 1994, ApJ, 437, 83 
Barnes, J. E. 1988, ApJ, 331, 699 

Barnes, J. E. & Hernquist, L. 1992, ARA&A, 30, 705 
Breitschwerdt, D., Voelk, H. J. & McKenzie, J. F. 1991, A&A, 245, 
79 

Breitschwerdt, D. & Schmutzler, T. 1999, A&A, 347, 650 
Bruzual, A. G., & Chariot, S. 1993, ApJ, 405, 538 



Buote, D. A. 2000, MNRAS, 311, 176 

Carr, B. J., Bond, J. R. & Arnett, W. D. 1984, ApJ, 277, 445 
Calzetti, D. & Heckman, T. M. 1999, ApJ, 519, 27 
Gen, R., Miralda-Escude, J., Ostriker, J. P. & Rauch, M. 1994, ApJ, 
427, L9 

Gen, R. & Ostriker, J. P. 1999, ApJ, 519, L109 

Ghiao, R. Y. & Wickramasinghe, N. G. 1972, MNRAS, 159, 361 

Gowie, L. L., Songaila, A., Kim, T., & Hu, E. M. 1995, AJ, 109, 1522 

Gowie, L. L. & Songaila, A. 1998, Nature, 394, 44 

Dave, R., Hellsten, U., Hernquist, L., Katz, N. and Weinberg, D. H. 

1998, ApJ, 509, 661 
Dave, R., Hernquist, L., Katz, N., & Weinberg, D. 1999a, ApJ, 511, 

521 



33 



Dave, R., Gardner, J. P., Hernquist, L., Katz, N., & Weinberg, D. 

1999b, to appear in The Evolution of Galaxies on Cosmological 

Timescalps, pHs .1 F Reckman and T. J. Mahoney, ASP books 

(1999); istro-ph/991022CH 
David, L. P., t'ornian, VV. & Jones, C. 1990, ApJ, 359, 29 
David, L. P., Forman, W. & Jones, C. 1991, ApJ, 380, 39 
Davies, J. I., Alton, P., Bianchi, S. & Trewhella, M. 1998, MNRAS, 

300, 1006 

Davis, D. S., Mulchaey, J. S. & Mushotzky, R. F. 1999, ApJ, 511, 34 
Dekel, A. & Silk, J. 1986, ApJ, 303, 39 

Duley, W. W., Jones, A. P. & Williams, D. A. 1989, MNRAS, 236, 
709 

Efstathiou, G. 2000, MNRAS, 317, 697 

Ellison, S. L., Songaila, A., Schaye, J. & Pettini, M. 2000, AJ, 120, 
1175 

Ferrara, A., Aiello, S., Ferrini, F. & Barsella, B. 1990, A&A, 240, 
259 

Ferrara, A., Pettini, M., & Shchekinov, Y. 2000, MNRAS, 319, 539 
Ferrara, A. & Tolstoy, E. 2000, MNRAS, 313, 291 
Finoguenov, A., David, L. P., & Ponman, T. J. 2000, ApJ, 544, 188 
Folkes, S. et al. 1999, MNRAS, 308, 459 

Fukugita, M., Hogan, C. J. & Peebles, P. J. E. 1998, ApJ, 503. 518 

Fukumoto, J. & Ikeuchi, S. 1996, PASJ, 48, 1 

Gavazzi, G., Pierini, D. & Boselli, A. 1996, A&A, 312, 397 

Gnedin, N. Y. & Ostriker, J. P. 1997, ApJ, 486, 581 

Gnedin, N. Y. 1998, MNRAS, 294, 407 

Gould, A. , Flynn, C. & Bahcall, J. N. 1998, ApJ, 503, 798 

Granato, G. L., Lacey, C. G., Silva, L., Bressan, A., Baugh, C. M., 

Cole, S., & Frenk, C. S. 2000, ApJ, 542, 710 
Gunn, J. E. & Gott, J. R. I. 1972, ApJ, 176, 1 
Haardt, F. & Madau, P. 1996, ApJ, 461, 20 
Haiman, Z. & Loeb, A. 1997, ApJ, 483,. 21 
Heckman, T.M. 2000, |astro-ph/000907E| 

Heckman, T. M., Rotiert, (J. , l^eitnerer, C. , Garnett, D. R. & van 

der Rydt, F. 1998, ApJ, 503, 646 
Heckman, T. M., Lehnert, M. D., Strickland, D. K., & Armus, L. 

2000, ApJS, 129, 493 [HLSA] 
Heckman, T. M., Armus, L. & Miley, G. K. 1990, ApJS, 74, 833 
Hellsten, U., Dave, R., Hernquist, L., Weinberg, D. H., & Katz, N. 

1997, ApJ, 487, 482 
Hellsten, U., Hernquist, L., Katz, N. & Weinberg, D. H. 1998, ApJ, 

499, 172 

Hernquist, L. 1992, ApJ, 400, 460 
Hernquist, L. 1993, ApJ, 409, 548 

Hernquist, L., Katz, N., Weinberg, D.H. & Miralda-Escude, J. 1996, 

ApJ, 457, L51 
Hernquist, L. & Katz, N. 1989, ApJS, 70, 419 
Impey, C, & Bothun, G. 1997, ARA&A, 35, 267. 
Jones, A. P., Tielens, A. G. G. M., HoUenbach, D. J. & McKee, C. 

F. 1994, ApJ, 433, 797 
Katz, N. ,' Weinberg^ D. H. & Hernquist, L. 1996, ApJS, 105, 19 
Kennicutt, R. C. 1998, ApJ, 498, 541 

Kim, S. -H. , Martin, P. G. & Hendry, P. D. 1994, ApJ, 422, 164 
Kobayashi, C. & Arimoto, N. 1999, ApJ, 527, 573 
Laor, A. & Draine, B. T. 1993, ApJ, 402, 441 
Lehnert, M. D. & Heckman, T. M. 1996, ApJ, 462, 651 
Lehnert, M. D., Heckman, T. M. & Weaver, K. A. 1999, ApJ, 523, 
575 

Mac Low, M. -M. & Ferrara, A. 1999, ApJ, 513, 142 | 

Madau, , P ., F errara, A., & Rees, M. 2000, ApJ, submitted; astro- 

ph/c |010158| 

Madau, P. & Pozzetti, L. 2000, MNRAS, 312, L9 
Martin, C. L. 1999, ApJ, 513, 156 
Mathis, J. S. & Whiffen, G. 1989, ApJ, 341, 808 
Metzler, C. A. & Evrard, A. E. 1994, ApJ, 437, 564 
Metzlet_-C — A — Evrard, A — E — 1997, ApJ, submitted 



ph/E 710324 

Murakami, 1. & Babul, A. 1999, MNRAS, 309, 161 



tstro. 



ApJ, submitted; a,stro-ph/010628: 



^eir 

i 



Dave, R. 2001, 



Mushotzky, R., Loewenstein, M., Arnaud, K. A., Tamura, T., 

Fukazawa, Y., Matsushita, K., Kikuchi, K. & Hatsukade, I. 1996, 

ApJ, 466, 686 
Nath, B. B. & Chiba, M. 1995, ApJ, 454, 604 
Nath, B. B. & Trentham, N. 1997, MNRAS, 291, 505 
Ostriker, J. P. & Gnedin, N. Y. 1996, ApJ, 472, L63 
Pagel, B. E. J. 1997, Nucleosynthesis and Chemical Evolution of 

Galaxies (Cambridge: Cambridge University press) 
Pagel, B. 1999, invited review, Ringberg Workshop: 'Galaxies in the 

Young TTnUrovaoTT' Aug 9-fi 1999, Hans Hippelein (ed.), Springer- 

Verlag; istro-ph/9911204 
Pearce, F. R., Jenkms, A., Frenk, C. S., Colberg, J. M., White, S. 

D. M., Thomas, P. A., Couchman, H. M. P., Peacock, J. A., & 

Efstathiou, G. 1999, ApJ, 521, L99 
Pei, Y. C, Fall, S. M. & Hauser, M. G. 1999, ApJ, 522, 604 
Persic, M. & Salucci, P. 1997, ASP Conf. Ser. 117: Dark and Visible 

Matter in Galaxies and Cosmological Implications. 
Pettini, M., Ellison, S. L., Steidel, C. C. & Bowen, D. V. 1999, ApJ, 

510, 576 

Pettini, M. , Steidel, C. C., Adelberger, K. L., Dickinson, M. & 

Giavalisco, M. 2000, ApJ, 528, 96 
Pettini, M. , King, D. L., Smith, L. J. & Hunstead, R. W. 1997, ApJ, 

478, 536 

Pettini, M., Shapley, A., Steidel, C. C, Cuby, J.-G., Dickinson, 
M., Moorwood, A.F.M., Adelberg er, K , L ,, Dickii jison, M., & 



Giavalisco, M. 2001, ApJ, in press; istro-ph/0102456 
Prochaska, J. X. & Wolfe, A. M. 2000, ApJ, 533, Lb — 
Quilis, v., Moore, B., & Bower, R. 2000, Science, 288, 1617 
Rauch, M., Haehnelt, M. G. & Steinmetz, M. 1997, ApJ, 481, 601 
Renzini, A. , Ciotti, L. , D'Ercole, A. & Pellegrini, S. 1993, ApJ, 419, 

52 

Renzini, A. 1997, ApJ, 488, 35 

Renzini, A. 1999, to appear in "Chemical Evolution from Zero to 
,Righ Redshift", ed. by J. Walsh & M. Rosa (Berlin: Springer); 



^stro-ph/990236l| 
Renzmi, A. ZWU, to appear in "Large Scale Structure in the X-Ray 
. Genrgantnpnulns & M. Plionis (Gyf sur Yvette, 



Universe", ed. I 



Dickinson, M. & 



Ed. Frontieres); istro-ph/0001312 
Salucci, P. & Persi c, M. 1999, MNR AS, 309, 923 
Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11. 1 
Scannapieco, E. & Broadhurst, T. 2001, ApJ, 549, 28 
Schaye, J., Rauch, M., Sargent, W. L. W., & Kim, T. 2000. ApJ, 
541, LI 

Shustov, B. M. & Vibe, D. Z. 1995, Astronomy Reports, 39, 578 
Simonsen, J. T. & Hannestad, S. 1999, A&A, 351. 1 
Songaila, A. 1997, ApJ, 490, LI 
Songaila, A. & Cowie, L. L. 1996, AJ, 112, 335 
Steidel, C. C, Adelberger, K. L., Giavalisco, M., 

Pettini, M. 1999, ApJ, 519, 1 
Steidel, C. C, Adelberger, K. L., Shapley, A. E., Pettini, M., 

Dickinson, M. & GiavaUsco, M. 2000, ApJ, 532, 170 
Strickland, D. K. & Stevens, L R. 2000, MNRAS, 314, 511 
Tegmark, M., Silk, J., & Evrard, A. 1993, ApJ, 417, 54 
Tegmark, M., Silk, J., Rees, M. J., Blanchard, A., Abel, T. & Palla, 

F. 1997, ApJ, 474, 1 
Theuns, T., Mo, H. J., & Schaye, J. 2001, MNRAS, 321, 450 
Wang, B. & Heckman, T. M. 1996, ApJ, 457, 645 
Weinberg, D.H., Dave, R., Gardner, J. P., Hernquist, L., & Katz, N. 
1999, in "Photometric Redshifts and High Redshift Galaxies", eds. 
R,. Weymann, L. Storrie-Lombardi, M. Sawicki & R,. Brunner, (San 
Francisco: ASP Conference Series), p. 341 
Weinberg, D.H., Hernquist, L., & Katz, N. 2000, ApJ, submitted; 

a.stro-ph/0005340| 
Wemgartner, J.U., & Draine, B.T. 1999, ApJ, 

ph/990725 
Wcingartnor, .J.C. 



submitted; astro- 



Draine, B.T. 1999, ApJ, submitted; astro- 



ph/0010117 
White, D. A. '2000, MNRAS, 312, 663 
Zaritsky, D. , Kennicutt, R. C. , Jr. & Huchra, J. P. 1994, ApJ, 420, 
87 

Zhang, Y., Anninos, P., & Norman, M.L. 1995, ApJ, 453, L57 
Zucca, E. et al. 1997, A&A, 326, 477 



